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SUMMARY 


The analysis of remote sensing data is very important in 
mapping surface and subsurface features. The knowledge of these 
features is helpful in estimating country’s hidden resourses. 
The LANDSAT MSS and TM data are available from National Remote 
Sensing Agency, Hyderabad. In the present work ,we have 
analysed MSS and TM data using image processing technique. The 
classified image has been compared with the toposheet prepared 
in the year 1976 by Survey of India. The various surface 
features such as river, roads, railway lines , soils, vegetated 
cover, forest areas, water bodies and built-up areas are 
classified using MSS and TM data. The information about a scene 
using MSS data is more because of larger area covered. However 
the resolution of the image obtained by TM data is higher 
compared to that of MSS imagery. Ue have taken remote sensing 
data of path number 144, row number 41 of MSS and quadrant ’D* 
of data and carried out digital analysis using Bayes technique 
and compared the results. He have found that both the data have 
their own limitations in mapping various surfacial features. 
The detailed analysis and superiority of these data has been 
discussed in the present work. 


CHAPTER I 


INTR<M>UCTI(»I 


1.1 Concept of Romofe Sensing 


The art of remote sensing started with the objective of 
locating the enemy camps and strategic locations. Remote 
sensing means getting the information about an object without 
coming in physical contact with it. This technique is based 
on the electromagnetic energy either reflected or emitted from 
the object. The wavelength range can be in visible, in the 
infrared or in the microwave range of electromagnetic 
spectrum. The components of a basic remote sensing system are 


shown in 

Fig. (1.1) 

. Conventional remote 

sensing 

uses 

the 

reflected 

energy 

in visible and 

infrared 

range 

in which 

sensors are used 

to detect the 

incoming . 

energy 

from 

the 

object. 

The technique has been 

used in 

a wide 

range 

of 

problems . 

Large 

data handling is 

made possible 

by 

the 


expanding computer network in different countries. The 
visual interpretation, assisted by aerial photos has been 
used as a complementary method. Amongst the satellite data 
available for the users, SPOT data look very promising giving 
EOffi resolution in the MSS mode and 10m in the panchromatic 
mode. Landsat-6 to be launched in near future is expected to 
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carry enhanced thematic mapper (ETM) . This shall incorporate 
additional band in the "panchromatic" band(0.5 to 0.86 m) 
having ground resolution of 15m. Sometimes the sensor is 
equipped to send the energy towards the object and 
getting the reflected energy. This type of sensing is known 
as active remote sensing and example is Radar. Radar 
sensors operate in the microwave frequency range (i.e.0.3 
to 300 GHz or 1 mm to 1m. range of wavelength) of 
electromagnetic spectrum. Microwave sensing also uses 
passive sensors like the microwave radiometers. Microwave 
sensors have the advantage of penetrating the atmosphere 
during day and night and also they can penetrate deeper in 
the ground. The common forms of microwave sensing are side 
looking airborne radar (SLAR), plan position indicator 
(PPl) and synthetic aperture radar (SAR). SLAR has been used 
for military reconnaissance. It has also been used for 
mapping the terrain in various countries, water resources and 
to prepare the geologic maps. Radar imagery has been used for 
the mapping of vegetation and crop identification in 
various countries. 

1.2 Remote SSensing Today 

Starting from 1960, remote sensing now is being used in 
almost all the branches of physical sciences. The major areas 
of application being geography, geology, civil engineering, 
forestry, meteorology, agriculture and oceanography. July 
1972 saw the launching of the first remote sensing satellite 
namely ERTS-1 by United States of America and subsequently 



the availability of data products in various countries 
assumed the future of the newly developing techniques. 
Since 1972, five such satellites have been sent for the 
purpose by United States of America , the latest being 
Landsat 5 (D*). SPOT, a remote sensing satellite of France 
became operational from 1985 giving resolution of as low as 
20 m. 

India has now joined the select band by sending IRS-1A 
(Indian Remote Sensing Satellite). The sensors work in the 
multi spectral scanner (MSS) mode to get the data in four 
bands. The details of 1RS-1A are given in table 1.1. 

1 .3 Remote Sensing Systems 

Figure (1.1) shows the basic elements of a remote sensing 
system. The main component of the system is the sensor which 
can be even a photographic camera for' the simplest of cases. 
At present the data is being collected by Landsat-5 of 
Landsat series. SPOT is also operational from 1985. The 
arbital characteristics of Landsat-4 and-5 are given below: 

[i) Sunsynchronus orbit. 

;ii) Pass at the equator of each orbit at 9.45 am. local sun 
time. 

iii) Near polar orbit, inclination of 98.2 with respect to 
the equator. 

iv) Altitude of 705 km. 


v) Repeat day period of 16 days. 
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Table 1.1 Indian Remote Sensing Satellite CIRS-1A> 


Objectives of Mission — To provide agricultural, geological and hydrologice| 

data for survey and management of nature 
resources 

Orbit Characteristics — Sunsynchronus, 904 km. altitude and having 22 . da 

repeat cycle. 

Sensor Applications Number of Spectral Resolution Swath 

Channels Range(pm) (m) (Km> 


LISS-I Landuse, urban 4 0.45-0.86 73 148 

(Linear Imaging planning, mapping. 

Self Scanning agriculture,water 

Sensor} resources, forestry, 

geological and 
mineral resources 

LISS-II -do- 4 0.45-0.86 36 148 



Energy 
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( 5 ) 

Real-time (instantaneous) 
data handling system 



( 6 ) 

Multiple data users 


Fig. 1 • i Cornponents o*f b , basic Reiniote Sensing Sys'ten'i 
(Lillesland and Kiefer 1983) 



Fig- 1-2 Sensors Onboard Landsat— 4 

(Lillesland and Kiefer 1985) 






The sa'telli'te collects the data by the following sensors— 

(i) Return Beam Vidicon 

(ii) Multi spectral scanner 
( iii IThematic mapper. 

Sensors onboard Landsat are shown in Fig. (1.S). 

The details of the MSS mode and thematic mapper are given 
Tables 1.S and 1.3 respectively. Details of SPOT satellite 

system are given in Table 1.4 

1-4 Remobe Sensing in Civil Engineering 

The application of remote sensing extends right from the 
terrain investigation to urban planning, in water resources, 
environmental studies, and a host of other aspects in which it 
acts as an aid to the more conventional methods. 


The remote sensing methods are highly cost effective 
means to get up-to-date and repetitive information^ Route 
mapping for highway location and design has been one of the 
most common projects in civil engineering. Application of 
remote sensing for site investigation emphasises on 


locating the features of geotechnical significance. 


Remote 



Table l.£ MSS Bands 


Landsat MSS Band No. Band Width Wave Band Name 

(^m) 

Landsat 1,2 & 3 Landsat 4 & 5 

5(D*) 

4 1 0.5-0. 6 Green 

5 2 0.6-0. 7 Red 

6 3 0.7— 0.8 Near Infrared 




Table 1.3 Thematic Mapper Bands 


Band 

No. 

Band Name 

Band Uidth 

(/jm) 

Characteristics 

1 

Blue/Green 

0.45-0.52 

Good water penetration 

strong vegetation 

absorbance 

2 

Green 

0.52-0.60 

Strong vegetation 

reflectance 

3 

Red 

0.63-0.69 

Very strong vegetation 

absorbance 

4 

Near 

Infrared 

0.76-0.90 

High land— water 

contrasts, very strong 

vegetation reflectance 

5 

Near-middle 

Infrared 

1.55-1.75 

Very mositure sensi- 
tive, effect of snow 

and clouds 

6 

Thermal 

Infrared 

10.4-12 

Very sensitive to soil 

moisture and vegetation 

7 

Middle 

Infrared 

2.08-2.35 

Good geological 

discrimination 


Table 1.4r SPOT Satellite System Specifications 


Orbital Details 


Altitude (km.) 

Period (min) 

Equatorial spacing of orbits (km.) 

Number of orbits/day 

Number of days for repeat coverage 

Equator crossing time 

Swath width (km) 

Sensor Details 

High Resolution Visible (HRV) Scanner- 
HSS mode — 

Number of spectral bands 
Spectral ranges (^Jm) 


Number of individual sensor detectors 
Ground Resolution 
Number of grey levels 
Date rate (mega bits/s) 

Black and white mode — 

Number of spectral bands 
Spectral range (^m) 

Number of individual sensor detectors 
Ground resolution 
Number of grey levels 
Date rate (mega bits/sec) 


832 

101 

1084 

14.2 

26 

10.30 

60 


3 

0.50-0.59 

0.61-0.68 

0.79-0.89 

3000 

20m 

256/(8 bit) 
25 


1 

0.51-0.73 

6000 

10m 

128 (6 bit) 
25 
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sensing plays important role in environmental engineering 
which includes its use in the study of the thermal and 
chemical pollution. Use of remote sensing in the studies of 
soil units is well established. Remote sensing in geology is 
used to give information on structure and lithology of rocks 
and features like f olding, faulting etc. Black and white, and 
colour infrared photographs have been used to map crops, crop 
growth etc. 

In the hydrological application surface and subsurface 
water sources can be inferred from the study of channels, 
lakes, drainage patterns etc. These features can be used detect 
ground water, which is indirectly inferred with the use of 
remote sensing. Other areas of study include flood and 
wetland studies. In areas of irrigation, large differences of 
soil moisture could be delineated to know the water use and 
stress on the crops. Areas of recent rainfall are detected 
in semi arid and arid regions. Another important area in 
which the remote sensing has been applied is in the 
distribution of snow fields and the likely volume of water 


from their melt 
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CHAPTER II 

FUNDAMENTALS OF ANALYSIS 

This chapter covers the properties of the cover surfaces 
and the reflectance response produced by different land cover 
types. A general description of pattern recognition 
techniques and the classifiers to extract useful information 
has been covered. Broad ideas of digital image processing and 
its usefulness in feature extraction is given in this 
chapter. 


2. 1 Reflectance from Cover Surfaces 

A graph of the spectral reflectance of an object as a 
function of wavelength is termed as spectral reflectance curve. 
The configuration of spectral reflectance curve gives insight 
into the spectral characteristics of an object and dictates 
the choice of the wavelength region(s) in which remote sensing 
data are required for particular application. The 

presence of peaks and lows in the reflectance curve of a 
particular type of feature characterise one surface from 
others. The spectral responses measured by remote sensing 
sensors permit an assessment of the type and/or condition of 
the surface features, and hence called spectral signatures. 
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Typical reflec'tance curves for vege'ta'tiorir soil and water 
are shown in Fig. (2.1> In case of water the most distinctive 
characteristic is the energy absorption in near infrared 
wavelengths. The reflectance properties of water bodies are 
modified by other factors such as sediments, bio-mass etc. The 
reflectance decreases with the increase in depth of water 
bodies in the near infrared wavelengths. Water containing 
large quantities of suspended sediments resulting from soil 
erosion normally have much higher visible reflectance than 
clear waters. The effect of increase in chlorophyll 
concentration is seen in blue wavelength where reflectance 
decreases. These features can be used to estimate the 
concentration of algae, presence of tanning dyes and the 
detection of a number of pollutants. 

Vegetation gives lower reflectance in range of 0.45 and 
0-67 /jm due to the pigment in the plant leaves. In the near 
infrared <0.70 to 0.8 /jm) the reflectance of healthy 
vegetation increases dramatically. Plant reflectance in the 
range of 0.70 to 1.3 ^im results from the internal structure of 
plant leaves. Since the structure is highly variable between 
plant species, reflectance measurements in this range help to 
discriminate between different types. Leaf reflection is 
inversely related to the total water present in a leaf. 
Reflectance is also controlled by the thickness of a leaf. 



Hemispherical reflectance (%) 
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Wavelength (fjim) 

Fig, 2.1 Typical Reflectance Curves for 
Water Soil and Vegitation 

(Paul J. Curran 1985) 



Fig. 2.2 Reflectance Curves for Soils 


(Paul J. Curran 1985) 



14 


Reflectance from soil is a function of particle 

size, texture, mineral composition, organic content, moisture 
content and temperature. Clay soils have strong structure in 
terms of packing and high moisture content, and hence 
give fairly low reflectance (presence of organic content 
also reduces the reflectance). Sandy soils on the other 
hand have a weak structure and present a fairly smooth 
surface and low moisture content and hence fairly high and 
often specular reflectance properties. In the visible 
wavelengths the presence of soil moisture considerably 
reduces the surface reflectance of soil. The same is true 
in near and middle infrared wavelengths also. Soil organic 
matter is dark and its presence decreases the reflectance 

from the soil upto organic content of 4-5X. For greater 

organic content further decrease in the reflectance is not 
- appreciable. 

One peculiar characteristic of soils is their 

association with vegetation. Often a particular type of soil 
is associated with a particular typ of vegetation. Also the 

extent of cover and the nature of vegetative cover (fully 

developed or less developed) affects the soil reflectance. 
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The effect of vegetation can be in the following ways— 

\ 

(a) due to shadow of vegetation 

(b) due to partial cover and 

(c) due to change in organic content and moisture content 
brought out by vegetative cover. 

The reflectance curves for soils are shown in Fig- (2.Z) . 


2.2 Pattern Recognition and Classification 

The remote sensing data is full of surface and subsurface 
information. This information is only extracted when analyse 
and interpret the data accurately. The analysis and 

interpretation is carried out by using pattern recognition and 
classification technique which is carried out either with 
visual interpretation or with computer. Visual interpretation 
involves the human experience and skill and the accuracy in 
the interpretation depends on the interpreter. The 

qualitative aspects of visual interpretation are - 
C15 Tone - In the satellite imagery there is variation of 
tone or shades, depending on the radiation reflected from 
different surfaces. It is expressed as light, medium and 
dark. 

Cll> Texture - This is defined as therate of change of tonal 
values, generally classified as coarse, medium and fine. 
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Cliii> Size - Depending on the scale of the imagery available, 
relative difference can be identified. 

Clv]> Shape - Different ground features have different shapes 
due to structure and topography which can be identified from 
the imagery. 

Cv3 Shadow- Shadows, cast particularly by the tall objects are 
used for the identification. 

CviJ Association— Some features on the imagery can not be 
directly identified as such but due to association with 
particular features their presence is indirectly inferred. 

CviiJ Pattern - Pattern is the spatial arrangement of various 
types of features. The arrangements can be in different 
forms. Geological patterns are either linear or curved 
depending on the geological features such as faults. Joints, 
folds etc. The topographical pattern is related to various 
types of drainage patterns (Fig. 2. 3). The drainage pattern 
and texture seen on any imagery are indicators of landform 
and type of bedrock and soil. The dendritic drainage 
pattern is a well integrated pattern formed by main stream 
with its tributaries branching and rebranching in all 
directions. The trellis pattern has one dominant direction 
and tributaries at right angles to it. The radial pattern is 
typical of volcanoes and domes. Anastomotic pattern with its 
meandering streams, cut-off meanders and ox— bow lakes is 
typical of the alluvial terrain. 



Lowlond orcos 




Hilly mounlaincous areas 


Freely developed Struclurallycontr- 




Anasiomotic 




Sub-dendrilic 



Sub-parallel jAnqular-anqulat^ 



Parallel Rccianqular 


Fig. 2.3 Drainage Patterns 
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2.3 CoinpurLer Aided Analysis and Digit/al Image Processing 

The interpretation of the remotely sensed data involves 
large data to be handled to get the useful information and 
hence the visual interpretation becomes difficult. Image 
classification using computer, automatically categorizes all 
the pixels in an image into number of classes which can be used 
to represent various themes of land cover, spatial pattern 
recognition the categorization of image pixels is done 
depending on the spatial relationship with surrounding pixels. 
The temporal pattern recognition of a scene can be used as 
an aid in feature identification. Image classification can 
be either supervised or unsupervised. 

Basic steps in the supervised classification are shown in 
Fig.<Z.4). In the training stage, sets of spectral data is 
collected to determine decision rules for the classification. 
In the classification stage each pixel]^ in the image data is 
categorized into the land cover class in which there is 
resemblance. After the categorization, the results are 
presented in the output stage. The unsupervised classifier 
does not utilize training data as the basis of classification 
and hence stage I is absent. 



Channel 


IMAGE DATASET 
(Five digital numbers^ 
per pixel) 




E: 

f- 


4- 

I 


p:i 

— 


III 

a— 


_ __ 


—r 

-_j__ 

H — 

— I 



— 

±: 

— 






CATEGORIZED SET 
(Digital numbers replaced 
by category types) 



(1) TRAINING STAGE (2) CLASSIFICATION STAGE (3) OUTPUT STAGE 


Collect numerical 
data from training 
areas on spectral 
response patterns 
of land cover 
categories 


Compare each unknown 
pixel to spectral patterns, 
assign to most similar 
category 


Present results 
maps 

tables of area data 
digital data files 


Fig. 2.4 Steps in Supervised Classification 
(Lillesland and Kiefer 1985) 



Fig. 2.5 Two Band Clustering of data 
(Lillesland and Kiefer 1985) 
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2. 4 Unsupervised Classllica'tlon 

The basic requirement for using such classifiers is 
that the different classes should be comparatively well 
separated. As a result of classification, analyst gets 
different spectral classes. These classified data have to be 

compared to some form of reference data to determine 

identity of spectral classes. The major advantage of such 
classifiers is that many of the classified classes might not 
be initially apparent to the analyst applying supervised 
classifier. When the spectral classes are too many it is 
advantageous to use such classifiers. The clustering 
algorithms may use 2 (Fig. 2. 5) or 4 channel data. K-class 
classifier is one such approach. In the present study, 

mode seeking algorithm has been used. 

2.5 Mode Seeking Algorithm 

Clusters are the homogeneous groups developed by the 
concentration of observation data around a point, called 
cluster centre. Any observation is assigned to the cluster 
only when the distance to the cluster centre is less than or 
equal to cluster threshold. In case when many cluster 
centres exist within the cluster threshold, the observation is 
assigned to the cluster with the closest centre. Due to the 
constant accretion in any cluster the cl ii^t er , the cluster 



centre is also modified. The step-wise procedure is as 
follows- 

(i) First data point is taken as the centre of first cluster. 

(ii) Calculate the square of the distance between the next 
data point and the first cluster centre. 

(iii) If the value so found is less than or equal to the 
cluster threshold then the point is assigned to the first 
cluster. 

(iv) If not, a new cluster with the above point originates. 

(v) With any data point the square of the distance of the 
data point with all the cluster centres is calculated and 
smallest of these values is chosen. If the selected distance 
is less than or equal to the cluster threshold the 
observation data is allotted to the nearest cluster. 

(vi) Grade which represents the members present in the cluster 
is constantly upgraded. 

The flow chart for the algorithm used, is given in Fig. (2.6). 

Stm& Supervised Classlflcat.lon 

To select the training data set for the supervised 
classifier complete knowledge of the geographic area is 
necessary. The quality and accuracy of the training data 


controls the 


accuracy 


of classification and hence 



H<k« th* first samplf tc b« tha 
cent'^r of th« first cluster 


Consider the next 
sample 


.ast l•mp1el 


Increase the cluster 
threshold 




Calculate the square of the 
distance from the sample to the 
centers of the existing clusters 
and hence find the nearest cluster 


^.^-^lusters ex-^ 
ceed the allowable 
'^'^0. of cl us- ^ 



^„,-''''^ampl e 
"^thln threshold of 
nearest ^ 
^V^cluster 


Assign the sample to the nearest 
cluster and update Its center as 
well as grade of membership 


Xreate • new cluster 
with the sample as 
th# center 




interpretation of image. The training data set has to be 
truly representative of the area. Sometimes two spectral 
classes might be required to adequately train a feature with 
large variation in spectral characteristics. In order to 
evaluate the training data, self classification of the 
training data is done to detect errors in the data. 


In the 

classification stage. 

actual 

classsif ication 

of 

the test 

sample is performed 

using 

mathemati cal 

and 


statistical appropaches. It is possible to use 2 or 4 channel 
data of MSS. Bayes decision rule is used for the 

classification purpose. This is an extension of maximum 

likelihood approach. The approach evaluates both the variance 
and covariance of the category spectral response patterns. The 
distribution for the clouds of points constituting a category 
is assumed Gaussian. With this assumption distribution 
of response pattern of any category can' be completely 
described by the mean vector and the covariance matrix. 
Statistical probability is calculated of a given pixel being a 
member of a particular category. The resulting probability 
density functions are shown in Fig.(S.7). Probability of each 
pixel is checked and the pixels are assigned to the most likely 
class. In the Bayes decision the probability of 

misclassifi cation is calculated for each category by assigning 
a “priori probability" for each class and assigning a weight 
associated with the error of misclassif ication. 




2. 7 Bayes Decision 


Any measurment X is assigned to the class G, with the 
highest conditional probability, minimizing the total error of 

• V 

classification. Amongst the G group or classes measurment X is 
assigned to group i if— 

P(G./X) > P(G^/X) for all 
where P(Gj^)= Probability of class "i" 
and P(Gj)= Probability of class "j". 

The final form of Bayes rule is 

P(X/G.) P(G.) > P(X/Gj) P<Gj) for all js'i 
To apply the Bayes rule we have to use the value of P(G^ ) and 
P(X/G^) for each group. Now assuming that within each group the 
varibles that make up the measurment vector X have a 
multivariate normal distribution, the form of condition^ji 
probability P(X/G. ) becomes 

P(X/G.)= 1 expC-1 (X-/J.)E (X-#j.)) 

where mean of class or group *i* 

variance covariance matrix of class 'i* 
and P(Gj)= A priori probability of class 'i*. 

Hence the estimation of P{X/G^) becomes calculating the 

parameters, the group mean vector and 2- the group 

covariance matrix. If there are n measurments in X (n bands) 

then there are 'n* group means in the mean vector and n(n+1) 

2 

elements in the covariance matrix making a total of n -t-Sn 


quantities to be estimated for each group, but this is 
considerably less than the number of estimations required for a 
direct application of the Bayes rule. 


Using the normal form of P(X/Gj) we have 


Assign X to G. if— 


P(G- ) 


-1 






i.y'S 


exp C-1 (X-Ai.)E (X->i-)'n > 

2 ^ i ^ 


P<G^) 




rixa 




exp C-1 (X-^.)E 

2 J j -» 


for all 

Taking the normal log(ln) of both sides 

T 

-1 n ln<En)-1 ln|E-l “1 <X->j.)E- <X-^.)'+ In P<G, > > 

" ■ ' g” ' g 1 — " g XIX 1 

T 

-1 n ln(2n)-1 InfE.I-l tX-*J.)Ei (X-jj.) + In P(G. ) 

" "■ 5 " ' ^ J g— J J J J 


for all ip'j 


Cancelling the common terms, 

-1 ... 


-ln|E- |-tX-#j. )E- (X-/J.) + In P(G, ) > 

X X 1 X 1 

” * x 

ln|Ej!-CX-Mj)Ej In PtG^) 

for all i^j 


■* ln|Eil+<X-AJ. )Ei <X-p.)^- In P(G.) < 

“*■ T 

InlE-l +<X-^i.)E. <X-TJ.> - In P(G.) 

Jl J J J J 

for all j^i 

T 

The quantity InjE^I +<X-# j ^ lEj ^X-jM^) “ In P(G. ) is referred to 

as "Discriminant Score” and 

d.(X) = lnlEil+<X-Mi>Ei <X-#J.r is called "Discriminant 


Function" 



As' a final rule of Bayes decision, assign the measurement *X 


to the group with the smallest value of dj<X>— In P(G^). 

Hence, assign 'X* to group i if— 

d.(X)-ln P(G.) < d.(X)-ln P<G.) for all j;-i. 

1 1 J J 

X is multivariate normally distributed in each of the 
classes . 

tiaximum likelihood Bayes classsifier is the most 
efficient classifier for categorization of the image pixels. 
For further study it is proposed that the urban area from the 
scene could be zoomed to get more detailed information. 

For the detailed vegetation analysis Kauth and 
Thomas(ref .4) have derived a linear transformation using four 
bands of landsat MSS. Crist and Cicone(ref.4 ) extended the 
same concept to Landsat TM data. For vegetation analysis and 
crop identtif ication three indices have been defined which 
are mathematically expressed as (Lillesand and Kiefer, 1985). 

(i) Vegetation Index (VI) = DNS - DN1 

_ DNS-DNi 

(ii) Normalized Vegetation Index (NVI) = t)N&+t>Ni ^ 


where DNS is reflectance value in near infrared band and DN1 is 
the reflectance value in visible band (i.e. 0.58 to 0.68Fin> 



1/s 


(iii) Transformed Vegetation 
Index (TVI) 

Here DN3 and DN4 are the reflectance values in band 3 and band 4 
of TM. 

Vegetated areas invariably give high value of TVI. 
Programme to compute the TVI for the whole scene to generated 
the TVI image is given (Index.fort) .Using this programme as 
many as six vegetation categories can be analysed. 

Z.6 Dlglbal Interactive Image Processing 

Digital image processing encompasses a wide variety of 
techniques and mathematical tools. A digital image is 

accepted as input and the image in some aspects is enhanced 
with image enhancement process. Image enhancement 

includes manipulating contrast, removing geometric 

distortion, edge smoothing and sharpening etc. A standard 
imagery is a general purpose product and cannot be 

controlled in a way users like hydrologists, geologists 
and urban planners desire. Digital images in the other hand 
can be made to fit for the specific purposes. To make the 
feature extraction easier standard false colour composites are 
prepared. The HSS false colour composite is generally created 
by exposing three of the four black and white bands through 


DN4-DN3 
DM4-^Di'ia ■ *‘’-= 
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different colour filters on to a colour film. In the false 
colour composite healthy vegetation appears bright red 
rather than green, water appears black, while 
sediment— charged water is powder blue. An urban area appear 
blue or grey blue. 

Image display is an integral part of image processing. 
The computer line printer is the simplest. The computer 
displays the image in the form of shade print on a line 
printer. One band density slicing (programme Den. FOR) 
can be used by coding the pixels with preassessed values and 
intervals. Different characters are used to represent 
the different density ranges. The advantage of this with 
method is due to its low cost. CRT is widely used and special 

OC 

high resolution CRT is specifically used for image pressing. 
The third type of device is the digital film recorder, which 
gives image of high quality. 

For image processing CAD-P facility was used. From 
the available grey levels of 0—255 in the data record, this 
can take up a maximum of 16 grey levels. Thus, the entire 
range of 0-255 has to be divided into 16 levels. The 
program developed by K.V. Rao(1988> is used in the image 
generation. A histogram analysis is tentatively carried out to 
slice the entire range to give good quality image. The system 
accepts in the maximum format of 480x640 pixels. The colour 
coding facility uses seven primary colours and a total of 
4096 colours in combination*. This is used to highlight 
specific features and zones. 





CHAPTER III 

DATA ACf^SITIOM AMD DCTAIDS OF AREA 

toltW 

This chapter deals ^ the procedure of data acquisition and 
details of the data products and their utility to carry out 
remote sensing analysis. 

3. 1 Details of Data 


To 

get 

the data 

for 

any area one needs 

to know the 

location 

of 

the area. 

A 

scene 

obtained by 

remote 

sensing 

satellite 

is 

located by 

certain 

path number 

and row 

number. 


Fig. (3.1) shows the Landsat pass over India and the path 

number, row number required for acquiring the data. Landsat 
data obtained for the present study consist of computer 

compatible tape for multi spectral mode (MSS) of iSth Sept. 
1986, which gave a cloud free picture. Thematic mapper 

data of quadrant *D* of the scene taken on 10th Nov. 1985. 

Paper ptTints of 1*250,000 in MSS for band S and band 4 and 
paper prints of 1*S50,000 in TM for band 
are also acquired. 


5, band 6, band 7 
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Fig. 3.1 


Landsat pass over India 



An area of 185 km x 185 km is covered in one Landsat 
scene. The information is obtained pixel by pixel and scan 
line by scan line. Multi spectral scanner scans the data in 
four bands and one pixel has the dimension of 79x79m. For 
thematic mapper the pixel has the dimension of 30mx30m except 
in band 6 (thermal Infrared) where the reolution is 120 m. The 
scanner scans the data in seven bands. 

The data collected by the scanning mirror are compared 
with the standard calibration source and sent to the detector 
where the intensity of reflectance is observed. Equivalent 
electrical signal ressults from the available reflectance, 
which is digitized and stored in computer compatible tapes 
(CCT), in binary form. The data of the Landsat— 4 MSS is in 
BIL (Band Interleave) format shown in Fig. (3.2). For the 
comparative study the thematic mapper data is of Landsat— 5 of 
loth Nov. 1985 for quadrant 'D* with path No. 144 and row 
No. 041. The quadrant identification for thematic mapper is 
shown in Fig. (3.3). 

Thematic mapper gathers information in seven bands and 
hence the volume of data becomes nearly eight times the MSS 
volume for one scene. The format for thematic mapper is also 
band interleaved (BIL); The data for each quadrant is covered 
by 2848 scan lines and 3056 pixels per line. One data product 
(i.e. of a quadrant) occupiefs 3 tapes in 1600 bpi density. 
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Hence the volume in terms of the tapes required becomes 1S 
times the MSS volume for one complete scene. The 3 volumes 
of each quadrant cover 1000, 1000 and 848 lines respectively. 
The format of superstructure is shown in Fig. (3.4 ). 

The data required for the study, i.e. CCT for MSS of 

Landsat— 4 and TM data for Landsat— 5 were acquired from NRSA 

B B 

B Hyderabad. Topo sheets 63B, 63-^ and 63 — ^ were acquired 
from the office of Survey of India, Lucknow. The 
computer compatible tape (CCT) for MSS originally acquired was 
in 32 bit. DEC-1090 can read the tape only in 36 bit 
hence the conversion is necessary. The conversion from 32 bit 

to 36 bit format was done CMC Ltd. Secunderabad. Volume 3 

A 

data of TM was used in ND— 560 system available in CAD— P 
for the generation of images. 

3.2 Detallfi ot the Area 

The area selected for the study is Lucknow area. The 
area enclosed by latitudes and longitudes is shown in Fig. (3.5). 


Physiography of" the Area 

The regional slope of the area is towards southeast and 
Ramganga, Gomti and Sarda rivers flow in the same direction. 


Tarai plains are well developed in the region. 


Gomti is a 
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PIG 3.4 tm format superstructure 



sluggish stream with intricate series of meanders- To the 


north of the area is river Ghaghra. Mean elevation of Lucknow 
is 111 m above mean sea level. . 

Weal, her and ClimaLe 

Typical indicators of climate are as under- 
(1) Mean annual ranifall about 100 cm. 

(E) (a) Mean temperature of the day - 18*C (January) 

(b) Mean temperature of the day - 30 C(April) 

(c) Mean temperature of the day — 30*C(July) 

(d) Mean maximum temperature - 40*C 

(3) Mean relative humidity 805t (at 8.30 hrs) in January. 

Mean monthly values of temperature and rainfall are shown 

against each month. 

Maximum temperature — E3.3 C 
January Minimum temperature ~ 8.4 C 

/■ 

Rainfall - 17.5 mm 

Maximum temperature — 25.9 C 

February Minimum temperature — 10.8 C 


Rainfall 


— 20.6 mm 



March 


April 


May 


June 


July 


August 


Maximum temperature 
Minimum temperature 
Rainfall 

Maximum temperature 
Minimum temperature 
Rainfall 

Maximum temperature 
Minimum temperature 
Rainfall 

Maximum temperature 
Minimum temperature 
Rainfall 

Maximum temperature 
Minimum temperature 

Rainfall 

Maximum temperature 
Minimum temperature 
Rainfall 


32. 7* C 
15.9*C 
8.6mm 

38.6* C 
21. 6“ C 
9. 1 mm 
40.8*C 
25.7* C 

% 

15.0 mm 
37.9*C 
27. 6* C 
88.9 mm 
33.6*C 
26.4*C 
308 . 1 mm 
32.5*C 
25.9*C 
286.5 mm 



Se ptember 


October 


November 


December 


Maximum temperature — 33.3 C 


Minimum temperature - 24.7 C 


Rainfall 


- 213.1 mm 


Maximum temperature — 33.0 C 

Minimum temperature — 19.2* C 

Rainfall — 35.1 mm 

Maximum temperature — 28.8*0 

Minimum temperature — 12.3* C 


Rainfall 


— 5.6 mm 


Maximum temperature — 24.4 C 


Minimum temperature — 8.5 C 


Rainfal 1 


6.3 mm 


Geology of bhe Area 

The detailed geology of area is not very clear. There is 
a thick bed of alluvium, mostly sand and silt. The depth of 
alluvium is of the order of 1000 m. Beneath the alluvium 
cover are the unconsolidated siwalik sediments# and older 
tertiary formations. Underlying these are the more 

consolidated older formations such as upper Gondwana and 
Cretaceous. In the porous beds there are large reservoirs of 
fresh water. Water table in the area is upto 100-^50 m 

below the ground level- Rock outcrop is absent in the area. 



Soil Typres 


Detailed soil investigation has been carried out in the 
area by soil survey and soil work division. Department of 
Agriculture U.P. Soil of the area is classsified in 8 
groups from the agricultural point of view. The different 
soil groups are— 

(1) Gomti Tarai 

(2) Gomti Flats 

(3) Gomti Flats (Halomorphic) 

(4> Gomti Uplands 

(5) Gomti Lowlands 

(8) Gomti Lowlands (Halomorphic) 

(7) Sai Uplands 

(S> Ravines 

From the engineering consideration the soil cover is all 


alluvium 



CHAPTER IV 


ANALYSIS AND INTERPRETATION 


Remotely sensed data product are available in the form of 
computer compatible tapes (CCT), imageries in different bands 
and negative films which are used for visual interpretation. 
Analysis with the help of computer compatible tapes is a cost 
effective means since it can be used to generate the images 
using image processing. Also different classification 
techniques can be attempted. For the analysis of data with the 
help of computer, data stored in CCT (Computer Compatible 
Tape) is used and the study area is selected. The preliminary 
stages are reading the tape and calculating the constants for 
the area. 

•4.1 Preliminary Analysis 

Stage I — Reading the Tape 

4 

In DEC-1090 system, it is only possible to read the tape 
in 36 bit format, original tape of the area processed in VAX 
at NRSA, Hyderabad had to be converted to 36 bit from 32 ^bit 
format by inserting dummy bytes. A programme to read the data 
is given (Record.FOR) . The actual data of the area start from 
the 3rd file. The first two files are Annotation and Header 
records respectively. These files have to be be skipped while 
reading.. To read the tape with the help of the programme, 
the details required include the density of the magnetic tape, 
the record length of each file. (The record length of first 
two files differs from the rest of the data records). 



The system available at CAD-P is ND-560 and it is 
possible to read the tape in 32 bit format. The tapes used 
are of 1600 bpi density. Although use of higher density tapes 
gives some problem (ie. more than 1600 bpi), it is possible to 
read the tapes of 800 bpi also. For reading the TM data, 
programme developed by K.V.Rao( 1988) (programme RTM.FORT) is 
used. 

Stage H — Calculation of Constants of the Scene - 

The corner coordinates (latitudes and longitudes) for the 
area are computed from the available imagery for MSS. The 
origin of reference axes is so chosen as to place the area in 
the first quadrant. The origin is E79*00*00" and N26‘00’00" 
as shown in Fig. (4.1). Corner points 1,2,3 and 4 respectively 
repressent line number and pixel number of (1,1), (1,3240), 
(2340,1) and (2340, 3240). 

The programme to convert the geographical coordinates of 
the corner points to conical orthomorphic coordinates is given 
in Rampal (1982). The results after conversion are shown in 
table(4.1) . 

To get the line number and pixel number of any point from 
the conical orthomorphic coordinates simple transformation is 


used: 




E79:00:00 
N 2 6 •.00:00 


FIG. 4.1 AREA UWDER STUDY WITH REFERENCE A^ES 



able 4.i 


Geographical and Conical Orthomorphic Coordinates of The 
Scene 


Corner Geographical Coordinate Conical Orthomorphic 

Coordinate 


Longitude Latitude Longitude<X) Latitude(Y) 


79"5S*5.54" 28*27*10.43" 85040.47 273383.82 


81*43*12.36" 28*10*52.14" 267045.64 245608.78 


79*26*7.68" 26*42*14.78" 43302.85 78443.95 


81*16*3.59" 26*26*22.62" 


225996.62 


50892.18 


X 


A1x + B1y + C1 


Y = AEx + BEy + CE 

where x = line number 
y = pixel number 

X = conical orthomorphic coordinate of longitude 

Y == conical orthomorphic coordinate of latitude 

A1, B1, C1, AE, BE, CS represent the six constants. 

To solve for the six unknown constants, at least six 
simultaneous equations are required. This requirement 
dictates a minimum of three ground control points needed. The 
three points chosen are the 3 corners of the imagery. 
The constants calculated are— 


A1 


-17.844S1 

B1 

s=: 

56.19178 

C1 

= 

8500S.1S 

AE 


-83.343E5 

BE 


-8.575189 

CS 

ss 

E73475.7 


The programme(Linpix.For ) uses these constants to 
calculate the line number and pixel number, from the given 
latitudes and longitudes of the points. 



A.2. VlsuAl Int«irprAt,ation rrom the Imagery 


The liSS imageries available for the study are taken by 
Landsat— 4 on i8th September 1986. The scale of the imageries 
is 1s250,000. For the comparison toposheets 63B, 63 B/13 and 
63 B/1 are used. 


Band 1 and band 2 of MSS are the best for studying the 
cultural centres like urban areas* road network etc. The 
foressted areas have high absorbance in band 2 and are easily 
identified. The vegetated areas show medium tone varying from 
one type to another. Sediment charged water can be detected 
in band 2. Band 3 and band 4 are good for delineating water 
bodies. 

River Gomti in the Lucknow area does not appear continuous 
in band 2 imagery due to lack of tonal contrast. Sarda canal 
going through the southern fringe of the Lucknow city 
appears in light tone. However* Haider canal is missed in 
the imagery. North eastern railway line has appeared in darker 
tone. Feature best identified on the band 2 imagery is of 
forest areas. Kukrail forest and other forest to the west of 
Lucknow city appear in very dark tone. Amongst the man made 
features road network is not identifiable in the imagery. 
Comparing it with band 4 imagery it is clear that the general 
information content is more in band 4. 



In band 4 imagery river Gomti appears in the dark tone 
and hence easily identified. The urban areas appear darker in 
tone than the surrounding features and hence easy to 
delineate. The urban area has grown considerably in the last 
few years. The expansion of the city has been some what 
limited in the southern part. There is considerable growth in 
the north and particularly the north-western part of Lucknow 
^rea. The road network is not easy to identify. Canals are 
clear in band E and the North eastern railway line connecting 
Kanpur is identified in medium tone due to the tonal contrast 
from the surrounding agricultural land. A recent bridge built 
over Gomti in the western Lucknow appears in the imagery. 

The MSS band 4 shows smooth texture , the variation in 
texture is very less. The tonal changes are apparent in the 
city area and the Kukrail forest area. The smooth texture 
and little variation in the surrounding forest areas seems 
to be agricultural land. The smooth texture can be due to 
alluvial deposits. Major geomorphic ' features are not seen 
in band 4 imagery. The soil appears well drained and medium 
textured in band 4. The presence of the same is apparent as 
point bar deposits near the banks in light tone. There is 
absence of shadow indicating the plain land. Scarcity of 
the drainage lines shows that the soil is quite permeable. 

The drainage pattern has intricate curves and meanders. The 

river has very stable course and is confined within the banks 
even in the greatest of discharges. In the imagery no 



abandoned drainage is identified. The tributary (Behta 
nadi ) Joining river Gomti appears well in the imagery. A 
segment of Kukrail nadi is clear in the imagery.. Amongst 
other features one cut-off channel and the Sarda canal (Lucknow 
Branch) have appeared. 


There is dramatic improvement in the thematic mapper 
imagery. The instantaneous field of view (IFOV) being 30mK30m. 
Thematic mapper imagery is available in band 5 taken on 7 
Sept. 1985 by Landsat-5. Some cloud patches are seen in band 5 
imagery. The effect is almost negligible in mapping surfacial 
features. On comparing with the MSS imagery it appears that 
the information content is more in TM imagery. The scale 
being same (i.e. 1:E50,000), TM band 5 imagery shows urban 
area in greater detail. Block like urban pattern is very 
clear. The new extension of the city in the north west 
direction has suburban features and is comparatively less 
dense. Major roads and road intersections are clear and have 
radial features. The Northern railway line which can not be 
seen in the MSS band 4, appears in TM band 5. The densest 
part of Lucknow, near Charbagh and the railway station are 
clear in TM imagery. All the other features indicated by 
band 2 and band 4 MSS are present. Unlike MSS band 4, some 
lines of secondary drainage are identified. This may be due 
to the greater detail coverage in TM. Also the time of 
satellite pass being September, an increase in discharge is 


expected with the saturation of soil. 



The band 6 TM imagery is not very useful from the point of 
view to study all the surface features of the area. The 
resolution being 120 m,many of the features nicely shown in 
band 5 are not clear, although major features can be 
mapped. The road pattern is not clear. The band 7 imagery of 
10th Nov. 1985, gives all the features mapped in band 5 but 
the effect of cloud is negligible. Road network and 
other features like newly constructed bridge are nicely 
shown. 


The improvement in the quality and information content in 

TM imagery is due to the fact that the resolution is greatly 

improved from 79 m to 30 m. Also the TM bands which are 

comparable to the MSS, have been more closely ranged 

(tre 

(table: 1.2) and theMre three additional bands included. 

A 

An analysis of band 2 in MSS of the area upstream of 
Kanpur, gives the features of river Ganga. The river has 
typical meandering anastomotic pattern (Fig. 2. 3). The river 
appears in medium light tone and the very light toned area 
represent the sediment deposits and silted water. The left 
bank is comparatively light toned than the right bank. This 
indicates that the area adj'oining the left bank is slightly 
higher land. The imagery is taken in September and there is 
evidence of parallel drainage due to large runoff. Shifting 
course of river is expected and one oxbow lake is 


identified. 



St,udy of Flood Plain 


There is no integrated drainage pattern due to level 
terrain and periodic flooding. The surface drainage is 
constituted by a principal stream and secondary streams in 
high runoff season. Exceptionally light toned areas are the 
point bar deposits iin the river bed. The oxbow lake in the 
darker tone suggests that it is occupied by organic soil and 
is little vegetated. Natural levees of the river areas are 
light toned. In the band 4 MSS imagery river Ganga and its 
local tributaries Isan nadi and Kalyani nadi are better 
marked. 


! !.T., 


RARY 


Acc. 


.No. A.iOjLib-tlvsU 



CHAPTER V 


RESULTS AMD DISCUSSIONS 

The remote sensing data were obtained in the -form of 
digital data and in the form of prints. The analysis of these 
data products has been carried out with the help of computer. 
The visual analysis of imageries have been also done. The 
computer aided analysis and visual analysis have been compared 
with the toposheets obtained from Survey of India. 

S. 1 Computer Aided Analysis 

The main advantage with the digital data available in the 
form of computer compatible tapes is that the analysis can be 
carried out in a desired way. The digital data can be used to 
classify the sample pixels in a given number of categories. 
Also, the data can be used to generate the digital images when 
requi red. 

In the present study two types of classifiers have been 
used. Mode seeking or clustering has been used to know the 
main spectrally separable classes in the area. This analysis 
is carried out for MSS data. The theory of clustering 
algorithm has been already explained in Chapter-II. 



If at a certain stage of a particular cluster. 

Number of samples in the cluster = n 
and cluster centre is <A,B) 

If a new point (C, D) is inserted in the cluster then the 
grade of the cluster becomes(n + 1) and the new cluster centre 
is also modified as 


f nA+C nB+D ^ 

n +1 ' n +1 * 

In the algorithm every time after the formation of a new 
cluster, the number of clusters are compared with the 
specified number of clusters. If the number of clusters 
available exceeds the specified number then the threshold is 
increased by a factor and iterations are performed. The 
procedure stops when all the data points have been classified 
in one or the other classes. All the possible two band 
combination have been tried. The results of the different 
combination of bands, i.e. band 1-2, band 1-3 etc. are shown 
in Table<5.1). 

It is to be noted that vegetation and built-up areas are 
not spectrally separable in each wavelength due to 
intermixing. Also the method is not efficient for the 
classification but gives the number of classes available in 
the data. Only band 1-3 combination gives some information 



Table 5.1 Results of Mode Seeking Analysis 



Number of 

Samples Correctly 

Classified as 

Band 

Combination 

Water 

Built-up 

Vegetation 

1-2 

50 


41 

1-3 

50 

49 

31 

1-4 

50 

50 

8 

2-3 

50 

48 

22 

2-4 

50 

50 

19 

3-4 

50 

50 

23 


where atleast three categories are clear, although good number 
of data samples from vegetation have been misclassif ied as 
built-up area. 


If we 

look 

at the results 

for 

different 

band 

combinations 

it is 

clear that water 

is 

classified 

very 


accurately in all the band combinations. Classification of 
built-up area is also coming good except in band 1— E 
combination. Vegetation is misclassif ied as built-up area. 
The misclassif ication is likely due to the recent built-up 
activity which is not represented in the toposheet of 1976. 
Another reason could be that large vegetated areas have been 
converted to built-up area. 


5.2 Results Using Bayes Decision 


This is a supervised parametric classification technique. 
This method is trained by a training data set which should be 
true representation of the category represented. To evaluate 
the decision, following parameters are calculated for 
training set — 

(1) Mean values of reflectance in various bands for 

different categories. 

(E) Variance of reflectance values in various bands for 


various classes 



(3) Covariance of reflectance values for all the combination 


of bands for all the classes. 

(4) Variance-covariance matrix for each classs. 

(5) Inverse of variance— covariance matrices. 

Test samples are tested by computing the probability 
and applying the decision rule. 

P(G^)= ^ Ei (x-p.)3 


-1 


In 


P(G. )= -1/2C(X-/li . ) 
1 1 


Ei 


(X-»i.)3= -1/2 


U 


For the conditional probability to be maximum, the value 
of *U* should be minimum. A given test sample is allotted to 
the class for which the value of ’U* is minimum.. 


5.3 Training the Data — 

For the training of data for Bayes classifier toposheets 
63B, 63B/1 representing Lucknow area are used. Training 
samples from river Gomti have been used to represent water 
body and likewise the training samples from the urban (built- 
up) area and Kukrail forest area etc. have been taken. The 
clustering analysis indicated that only three major categories 
could be successfully distinguished using MSS, hence training 



set (Table 5.2) comprises of four classes. Although, taking 
the training data in a grid pattern is better, it could not 
be done in such fashion due to narrow streak of Gomti and 
only localized forest area. The training set for TM (Table 
5.3 ) comprises of data in six classes. 

The test data set comprise of a total of 160 points in 
MSS, with 40 points each in one "predicted" class. For test 
sample of TM, 180 points have been taken. 

Firstly the training data itself is classified for MSS 
and TM. The results have come out with 100 percent accuracy. 
Tables (5.4) and (5.5) respectively show the confusion 

matrices for the training set data. Next, the test sample data 
is classified and the results obtained are presented in Tables 
(5.6) and ( 5.7 ). It is clear from the confusion matrix 

that the classifier has effectively classified the sample data 
for MSS and TM. This gives very accurate results for 4 
classes in MSS and 6 classes using TM data. For the MSS 

data water and urban area have been classified with 100 
percent accuracy. One sample each of the Forest/Vegetation 
class is falling under urban area and uncultivated land 

due to the similarity in reflectance response and absence 
of sharp boundary. In TM data water is classified with 100 
percent accuracy. For urban and sub-urban classes the accuracy 
obtained is 93.35t. Two samples from the urban area have 
been misclassif ied as sub— urban and vice-versa. This is due 
to the lack of sharp boundary between urban and sub-urban 
zones and a small localized pocket of denser sub-urban 


region. 



Table 5.2 Training Set for MSS 


Cat egor y-Wat er 


SI. No. 

Band 1 

Band 2 

Band 3 

Band 4 

1 

59 

53 _ 

52 

32 

2 

60 

52 

55 

31 

3 

62 

59 

60 

33 

4 

61 

53 

67 

41 

5 

59 

56 

65 

41 

6 

58 

54 

53 

36 

7 

60 

55 

56 

34 

8 

61 

50 

51 

26 

9 

62 

53 

65 

32 

10 

57 

58 

61 

35 


Category-Urban 


SI .No. 

Band 1 

Band 2 

Band 3 

Band 4 

1 

64 

66 

77 

51 

2 

67 

67 

80 

52 

3 

64 

57 

81 

55 

4 

67 

59 

69 

48 

5 

64 

62 

79 

52 

6 

64 

73 

75 

50 

7 

65 

67 

81 

61 

8 

66 

63 

78 

54 

9 

62 

58 

77 

52 

10 

65 

66 

79 

53 



Cat eg or y- Veget at 1 on 


SI. No. 

Band 1 

Band 2 

Band 3 

Band 4 

1 

48 

38 

99 

80 

2 

46 

36 

107 

87 

3 

48 

34 

107 

90 

4 

45 

37 

105 

85 

5 

45 

35 

106 

86 

6 

46 

37 

102 

87 

7 

49 

38 

102 

83 

8 

47 

36 

106 

84 

9 

46 

39 

105 

86 

10 

45 

32 

100 

89 


Category-Uncul tl vated 

SI. No. 

Band 1 

Band 2 

Band 3 

Band 4 

1 

62 

48 

126 

104 

2 

56 

52 

123 

102 

3 

58 

51 

130 

105 

4 

59 

52 

129 

103 

5 

61 

55 

131 

106 

6 

57 

52 

125 

105 

7 

60 

53 

124 

107 

8 

61 

54 

126 

107 

9 

58 

50 

132 

103 

10 

59 

56 

130 

101 


MSS TEST SAMPLE DATA 


S . NO . 

1 

2 

3 

4 

5 

6 

7 

8 
9 

•10 

•11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 


BAND1 

BAND2 

58 

53 

56 

53 

57 

58 

53 

52 

56 

57 

50 

55 

57 

57 

56 

54 

56 

53 

60 

57 

60 

53 

57 

52 

64 

57 

56 

56 

60 

46 

56 

56 

56 

57 

56 

53 

56 

65 

61 

62 

64 

57 

56 

57 

59 

57 

62 

55 

61 

53 

63 

57 

64 

56 

64 

57 

54 

50 

67 

57 

64 

59 

64 

57 

53 

52 

64 

60 

64 

54 

65 

54 

62 

55 

63 

54 

62 

55 

63 

54 

60 

57 

65 

60 

65 

62 

64 

57 

60 

57 

62 

58 

69 

69 

60 

72 

59 

62 

70 

65 

68 

65 

67 

68 

60 

68 


BAND3 

BAND4 

66 

42 

58 

38 

59 

34 

65 

46 

60 

38 

55 

36 

60 

38 

59 

36 

66 

41 

66 

42 

66 

42 

70 

46 

55 

29 

60 

32 

66 

33 

66 

35 

66 

38 

60 

33 

70 

47 

67 

36 

62 

33 

68 

42 

51 

33 

55 

33 

65 

41 

60 

34 

59 

35 

51 

29 

58 

33 

47 

23 

53 

33 

55 

29 

50 _ 

29 

61 ^ 

29 

55 

26 

60 

36 

61 

33 

63 

32 

60 

35 

61 

33 

81 

57 

82 

60 

80 

57 

82 

52 

83 

58 

79 

60 

80 

57 

81 

57 

78 

58 

79 

52 

84 

57 

90 

58 

88 

58 



54 

61 

55 

63 

56 

63 

57 

61 

58 

66 

59 

71 

60 

73 

61 

69 

62 

65 

63 

65 

64 

68 

65 

65 

66 

70 

67 

69 

68 

70 

69 

60 

70 

62 

71 

63 

72 

66 

73 

65 

74 

66 

75 

61 

76 

66 

77 

63 

78 

65 

79 

62 

80 

65 

81 

45 

82 

37 

83 

45 

84 

45 

85 

45 

86 

42 

87 

41 

88 

41 

89 

47 

90 

45 

91 

43 

92 

45 

93 

48 

94 

45 

95 

47 

96 

53 

97 

41 

98 

44 

99 

48 

100 

53 

101 

43 

102 

41 

103 

47 

104 

50 

105 

43 

106 

45 

107 

47 

108 

48 

109 

47 

110 

51 


70 

81 

63 

77 

64 

79 

70 

81 

64 

79 

69 

72 

70 

91 

59 

79 

65 

79 

66 

78 

68 

81 

58 

82 

69 

88 

72 

90 

67 

79 

55 

75 

56 

73 

65 

77 

65 

78 

55 

80 

60 

69 

62 

78 

68 

83 

65 

73 

69 

91 

72 

90 

67 

79 

34 

104 

46 

104 

32 

103 

44 

105 

39 

100 

42 

110 

28 

107 

53 

120 

48 

107 

32 

111 

37 

100 

36 

98 

81 

87 

29 

86 

37 

112 

49 

99 

28 

107 

45 

103 

57 

109 

49 

118 

37 

104 

28 

103 

42 

117 

40 

108 

44 

113 

50 

108 

37 

104 

41 

104 

49 

99 

48 

104 


56 

49 

59 

60 

55 

51 
61 

56 
53 

52 

53 
56 

59 

60 
58 

50 
50 

52 
50 

54 
49 
61 
65 
60 
63 
62 

53 
82 
78 
85 
81 
75 
87 

89 
84 

84 

90 
81 
82 
78 
81 

85 
81 
84 

84 
78 

85 
85 
84 

84 
92 

85 
81 

77 

78 
68 
78 



6 


111 

47 

1-12 

48 

113 

50 

114 

48 

115 

51 

116 

43 

117 

49 

118 

59 

119 

58 

120 

46 

121 

63 

122 

61 

123 

57 

124 

60 

125 

55 

126 

59 

127 

62 

128 

62 

129 

50 

130 

52 

131 

60 

132 

59 

133 

62 

134 

61 

135 

56 

136 

50 

137 

51 

138 

58 

139 

59 

140 

62 

141 

63 

142 

58 

143 

68 

144 

63 

145 

58 

146 

59 

147 

60 

148 

58 

149 

51 

150 

57 

151 

62 

152 

59 

153 

57 

154 

68 

155 

67 

156 

58 

157 

62 

158 

61 

159 

62 

160 

59 


45 

98 

69 

43 

98 

86 

40 

102 

82 

38 

109 

89 

36 

109 

92 

45 

79 

89 

32 

98 

80 

56 

98 

82 

58 

99 

100 

30 

112 

95 

49 

126 

105 

47 

125 

103 

53 

125 

102 

61 

126 

107 

54 

122 

106 

52 

129 

104 

62 

133 

108 

59 

129 

101 

55 

121 

100 

53 

119 

99 

57 

131 

102 

51 

113 

104 

55 

112 

108 

54 

125 

108 

56 

123 

100 

48 

108 

95 

49 

109 

92 

52 

129 

104 

53 

132 

103 

54 

131 

105 

54 

130 

106 

68 

110 

104 

58 

98 

103 

65 

133 

102 

54 

131 

107 

55 

130 

101 

56 

131 

105 

52 

120 

102 

53 

100 

103 

67 

129 

103 

63 

122 

111 

69 

128 

109 

47 

127 

108 

42 

122 

109 

47 

124 

106 

54 

128 

100 

68 

98 

89 

51 

127 

105 

52 

132 

103 

60 

130 

102 



63 


Table 5.3 Training Set for TM 


Category-Water 


SI. No./ 

1 

2 

3 

4 

5 

6 

7 

TVI 

bands 









1 

89 

41 

41 

28 

15 

124 

11 

55 

2 

93 

43 

43 

28 

20 

123 

8 

53 

3 

92 

40 

42 

33 

17 

125 

9 

61 

4 

93 

40 

42 

26 

18 

122 

10 

51 

5 

89 

40 

41 

25 

17 

121 

10 

50 

6 

93 

43 

45 

35 

18 

124 

9 

61 

7 

92 

41 

41 

25 

17 

124 

8 

50 

8 

89 

40 

41 

23 

17 

125 

9 

46 

9 

88 

40 

40 

24 

17 

121 

10 

50 

10 

89 

38 

39 

22 

21 

124 

10 

47 

Category-Urban 




• 




S.No./ 

1 

2 

3 

4 

5 

6 

7 

TVI 

bands 









1 

89 

38 

38 

34 

54 

127 

34 

66 

2 

88 

37 

38 

34 

52 

127 

33 

66 

3 

89 

37 

38 

38 

56 

127 

35 

70 

4 

90 

38 

38 

33 

53 

126 

35 

65 

3 

89 

38 

42 

36 

50 

127 

36 

65 

6 

89 

35 

36 

33 

45 

125 

34 

67 

7 

81 

35 

34 

38 

44 

124 

32 

74 

8 

84 

36 

36 

34 

46 

126 

33 

68 

9 

93 

40 

42 

35 

57 

126 

33 

64 

10 

86 

38 

38 

34 

52 

126 

34 

66 





Category-Sub -urban 


SI. No./ 

bands 

1 

2 

3 

4 

5 

6 

7 

TVI 

1 

86 

35 

35 

40 

59 

126 

40 

75 

2 

86 

36 

35 

44 

59 

124 

34 

78 

3 

86 

36 

34 

42 

60 

125 

30 

77 

4 

88 

36 

36 

40 

62 

127 

32 

74 

5 

88 

36 

38 

39 

59 

125 

31 

71 

6 

92 

41 

41 

38 

54 

123 

29 

68 

7 

91 

43 

46 

42 

61 

127 

37 

67 

8 

87 

37 

40 

41 

62 

124 

35 

71 

9 

86 

36 

36 

43 

58 

125 

34 

76 

10 

84 

37 

35 

42 

59 

126 

33 

76 


Category-Foresl 

SI. No./ 1 2 

bands 

3 

4 

5 

6 

7 

TVI 


1 

74 

30 

26 

48 

54 

122 

16 

89 

2 

75 

31 

29 

58 

45 

121 

16 

91 

3 

75 

30 

26 

59 

45 

121 

14 

94 

4 

78 

34 

32 

48 

55 

124 

15 

83 

5 

81 

33 

29 

50 

49 

127 

14 

87 

6 

75 

29 

26 

53 

47 

122 

16 

91 

7 

81 

36 

36 

50 

68 

126 

17 

81 

8 

76 

31 

27 

51 

52 

124 

16 

89 

9 

75 

30 

27 

52 

50 

127 

15 

90 

10 

75 

29 

26 

53 

47 

122 

15 

91 


Cat egory-Vegetat 1 on 


SI. No./ 

bands 

1 

2 

3 

4 

5 

6 

7 

TVI 

1 

89 

38 

38 

43 

63 

125 

35 

74 

2 

91 

41 

43 

45 

77 

125 

40 

72 

3 

90 

40 

40 

41 

66 

126 

41 

71 

4 

88 

38 

38 

62 

75 

126 

41 

86 

5 

93 

43 

49 

56 

93 

127 

36 

74 

6 

93 

44 

52 

58 

89 

127 

42 

74 

7 

90 

41 

42 

58 

77 

127 

40 

81 

8 

89 

41 

43 

50 

76 

125 

41 

75 

9 

86 

34 

35 

41 

66 

126 

41 

76 

10 

91 

41 

44 

51 

86 

127 

40 

75 


Category-Uncultivated 

SI. No./ 

bands 

1 

2 

3 

4 

5 

6 

7 

TVI 

1 

108 

54 

60 

56 

106 

126 

107 

68 

2 

110 

51 

57 

52 

103 

124 

108 

67 

3 

109 

53 

60 

55 

104 

122 

110 

67 

4 

97 

46 

52 

54 

102 

122 

112 

70 

5 

95 

46 

57 

56 

106 

126 

108 

70 

6 

99 

51 

52 

49 

100 

126 

109 

68 

7 

98 

49 

57 

58 

116 

125 

110 

71 

8 

102 

50 

58 

55 

99 

126 

111 

68 

9 

98 

51 

64 

67 

100 

120 

118 

72 

10 

97 

49 

60 

63 

114 

129 

115 

70 



54 

38 

42 

36 

53 

55 

35 

36 

34 

54 

56 

37 

39 

34 

56 

57 

37 

39 

33 

55 

58 

37 

37 

38 

49 

59 

35 

36 

36 

50 

60 

35 

35 

28 

38 

61 

37 

37 

38 

58 

62 

37 

37 

45 

74 

63 

38 

37 

40 

59 

64 

39 

39 

40 

60 

65 

37 

37 

38 

61 

66 

36 

34 

41 

59 

67 

41 

41 

39 

55 

68 

37 

37 

38 

60 

69 

36 

40 

40 

61 

70 

36 

36 

40 

59 

71 

37 

36 

41 

59 

72 

38 

36 

42 

60 

73 

38 

35 

38 

59 

74 

36 

36 

41 

43 

75 

41 

45 

42 

54 

76 

38 

39 

41 

57 

77 

38 

38 

39 

61 

78 

37 

38 

40 

61 

79 

41 

42 

41 

62 

80 

36 

36 

43 

59 

81 

40 

40 

40 

59 

82 

36 

37 

44 

59 

83 

37 

35 

41 

59 

84 

40 

41 

40 

60 

85 

36 

37 

43 

58 

86 

38 

43 

39 

60 

87 

37 

40 

40 

61 

88 

42 

44 

40 

68 

89 

40 

42 

42 

61 

90 

41 

44 

39 

58 

91 

35 

36 

37 

53 

92 

30 

26 

59 

44 

93 

30 

26 

52 

46 

94 

34 

32 

50 

46 

95 

29 

26 

49 

41 

96 

30 

29 

53 

50 

97 

35 

36 

53 

74 

98 

33 

31 

58 

64 

99 

32 

29 

53 

47 

100 

34 

32 

60 

59 

101 

36 

34 

54 

67 

102 

34 

34 

55 

51 

103 

31 

28 

53 

45 

104 

34 

35 

56 

53 

105 

30 

27 

51 

43 

106 

32 

27 

51 

47 

107 

30 

29 

58 

62 

108 

31 

27 

54 

45 

109 

33 

33 

48 

55 

110 

29 

25 

48 

4c 



61 


111 

32 

32 

49 

112 

32 

32 

50 

113 

30 

27 

58 

114 

34 

34 

58 

115 

36 

36 

50 

116 

35 

35 

50 

117 

34 

33 

42 

118 

34 

32 

57 

119 

33 

30 

53 

120 

35 

37 

44 

121 

37 

37 

50 

122 

42 

45 

49 

123 

40 

45 

45 

124 

40 

43 

46 

125 

43 

47 

50 

126 

40 

42 

47 

127 

30 

42 

49 

128 

44 

50 

62 

129 

41 

45 

54 

130 

40 

45 

47 

131 ' 

41 

43 

48 

132 

38 

38 

51 

133 

41 

44 

50 

134 

39 

41 

58 

135 

39 

40 

55 

136 

38 

40 

52 

137 

40 

45 

60 

138 

41 

45 

49 

139 

34 

34 

45 

140 

37 

38 

55 

141 

37 

37 

56 

142 

38 

41 

61 

143 

45 

52 

66 

144 

46 

55 

67 

145 

42 

45 

61 

146 

56 

65 

57 

147 

36 

36 

49 

148 

37 

38 

50 

149 

44 

51 

58 

150 

37 

37 

45 

151 

42 

43 

25 

152 

53 

60 

56 

153 

54 

58 

30 

154 

54 

61 

55 

155 

45 

55 

58 

156 

46 

59 

61 

157 

51 

68 

66 

158 

45 

55 

50 

159 

46 

53 

59 

160 

50 

60 

61 

161 

46 

45 

56 

162 

46 

55 

55 

163 

46 

51 

48 

164 

52 

61 

60 

165 

47 

58 

56 

166 

55 

66 

61 

167 

49 

61 

64 


46 

37 

44 

66 

62 

55 

36 

53 

48 

65 

63 

83 

77 
73 

78 
81 
68 
78 
85 
58 

65 

72 

70 

62 

60 

91 

77 

55 

85 
80 

86 
103 

91 

88 

120 

55 

45 

80 

70 

80 

100 

84 

102 

114 

114 

124 

94 

103 
106 

94 

84 

83 

99 

101 

94 

104 



69 


168 

46 

61 

57 

101 

169 

49 

57 

58 

106 

170 

49 

60 

58 

110 

171 

54 

67 

67 

123 

171 

51 

66 

60 

106 

172 

48 

59 

55 

100 

173 

49 

55 

58 

93 

174 

50 

60 

60 

99 

175 

45 

52 

50 

93 

176 

49 

59 

59 

89 

177 

47 

56 

58 

112 

178 

48 

57 

59 

100 

179 

47 

58 

54 

92 

180 

43 

49 

45 

69 
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Table 5.4 Confusion Matrix- MSS Training Data- 



Water 

Urban 

Forest/ 

Veg. 

Uncultivated 

Water 

10 




Urban 

— 

10 

— 

— 

Forest/Veg. . 

— 

— 

10 

— 

Uncultivated 

— 

— 

— 

10 


Table 5. 5 Confusion Matrlx-TM Training Data 

Water Urban Sub-urban Forest Veg. Uncultivated 

Water 10 — — — — 

Urban — 10 — — — 

Sub-urban — 10 — — 

Forest — — — 10 — 

Vegetation — — — — 10 

Unculti- 
vated 


10 




RESULTS FOR MSS 
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.NO. 

CLASS NO 

WU 

BU 

VU 

FU 

1 

1 

5.78 

39.80 

284.77 

2248.32 

2 

1 

7.23 

82.55 

364.86 

2583.72 

3 

1 

5.86 

71.10 

396.37 

2849.06 

4 

1 

33.77 

75.19 

230.64 

1831 .45 

5 

1 

5.99 

74.28 

357.15 

2563.88 

6 

1 

45.35 

172.60 

430 . 77 

2604.87 

7 

1 

3.42 

65 « S6 

356.83 

2594 . 75 

8 

1 

19.38 

27.32 

367.73 

2546.80 

9 

1 

11.42 

55.53 

288.93 

2240 . 42 

10 

1 

4.31 

26.84 

305.38 

2336.56 

11 

1 

5.22 

28.14 

293.32 

2309.35 

12 

1 

14.07 

38.46 

232.83 

1933.37 

13 

1 

12.76 

53.26 

491 . 17 

3503.13 

14 

1 

13.72 

80.45 

408 . 55 

2889.05 

15 

1 

20.36 

53.87 

442.90 

2761.62 

16 

1 

16.05 

89.75 

549 . 70 

2146.06 

17 

1 

10.54 

58.18 

317.14 

2422.55 

18 

1 

20.99 

104.80 

415.96 

2826.45 

19 

1 

13.50 

78.77 

404.41 

2810.05 

20 

1 

18.54 

48.45 

336.37 

1959.50 

21 

1 

10.17 

31 .71 

359 . 84 

2711.67 

22 

1 

9.68 

34.30 

411.68 

3071.00 

23 

1 

8.97 

50.18 

278.53 

2173.08 

24 

1 

4.72 

88.37 

470.22 

3177.12 

25 

1 

50.05 

88 . 67 

597.28 

3536.23 

26 

1 

5.62 

26.12 

311 .54 

2415.96 

27 

1 

15.09 

23.54 

351.91 

2499 . 26 

28 

1 

12.05 

33.19 

371.60 

2527.35 

29 

1 

17.65 

63.91 

522.52, 

3615. 12 

30 

1 

29.37 

105.04 

437.13 

2777.81 

31 

1 

38.72 

84.54 

634 . 69 

4264 . 06 

32 

1 

25.53 

52.81 

486.68 

3327 . 80 

33 

1 

12.76 

53.26 

491 . 17 

3503. 13 

34 

1 

30.47 

154.76 

537.52 

3198.03 

35 

1 

16.29 

48.01 

454.67 

3369.21 

36 

1 

8.54 

63.74 

532.34 

3671.15 

37 

1 

20.97 

32.75 

532.34 

3369.21 

38 

1 

1 .94 

40.25 

404.55 

3006.56 

39 

1 

3.73 

40.56 

413.38 

3045.44 

40 

1 

3.17 

37.59 

392.95 

2195.23 

41 

2 

31 .45 

11.72 

204.73 

1355.83 

42 

2 

79.79 

5.28 

324.68 

1367.42 

43 

2 

70.81 

1.76 

348.55 

1539.14 

44 

2 

22.20 

8.01 

213.95 

1595.35 

45 

2 

33.48 

12.04 

196.62 

1283.03 

46 

2 

58.83 

12.25 

269.73 

1329.69 

47 

2 

149.98 

9.87 

574.41 

1723.39 

48 

2 

62 . 62 

15.84 

492.44 

1489.39 

49 

2 

36.96 

21.80 

293.67 

1381.71 

50 

2 

119.62 

14.62 

462.70 

1960.57 

51 

2 

103.91 

10.20 

429.77 

1579.09 

52 

2 

92.70 

21.27 

443.10 

1434.25 

53 

2 

42.24 

16.11 

347.71 

1288.21 

54 

2 

56.15 

9.25 

440.45 

1531.91 

55 

2 

29.28 

3.21 

312.15 

1899.60 


DECISION 

•1 

1 

1 

1 

•1 

1 

1 

1 

1 

•1 

1 

1 

1 

1 

1 

1 

1 

1 

•1 

1 

•1 

1 

1 

1 

1 

1 

1 

•1 

1 

■1 

1 

1 

■1 

1 

1 

1 

1 

1 

1 

1 
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56 

2 

69.63 

7.17 

382.11 

1439.46 

2 

57 

2 

72.86 

14.73 

490.88 

1373.77 

2 

58 

2 

78.58 

1 . 12 

392.39 

1688.64 

2 

59 

2 

178.06 

16.97 

62 1 . 77 

2232 . 55 

2 

60 

2 

212.95 

57.09 

666.09 

1494.89 

2 

61 

2 

111.83 

11 . 14 

375.74 

1696.38 

2 

62 

2 

61 .38 

0.37 

378.28 

1756.13 

2 

63 

2 

62.76 

0.42 

397.09 

1830.18 

2 

64 

2 

102.79 

8.22 

474.43 

1839.67 

2 

65 

2 

55.80 

3.67 

271.25 

1515.03 

2 

66 

2 

150.14 

27.19 

556.52 

1524.45 

2 

67 

2 

149.80 

27.32 

611.06 

1452.97 

2 

68 

2 

168.78 

14.54 

571.19 

1710.04 

2 

69 

2 

15.57 

13.12 

217.07 

1747.88 

2 

70 

2 

23.67 

7.77 

253.93 

1850.31 

2 

71 

2 

44.30 

1 .80 

358.28 

1781 .20 

2 

72 

2 

60.96 

2.78 

384.66 

1943.66 

2 

73 

2 

49.41 

4.42 

247. 17 

1617.76 

2 

74 

2 

64.72 

7.20 

374.62 

2135.72 

2 

75 

2 

63.02 

20.40 

343.40 

1314.27 

2 

76 

2 

145.72 

15.83 

578.70 

1255.82 

2 

77 

2 

98.40 

27.44 

478.73 

1521 .87 

2 

78 

2 

97.49 

15.49 

483.48 

1183.87 

2 

79 

2 

79.44 

15.56 

507.24 

1195.86 

2 

80 

2 

68.58 

0 . 66 

416.76 

1775.79 

2 

81 

3 

387.81 

237.18 

17.55 

250.36 

3 

82 

3 

415.07 

356.72 

89.02 

285 . 78 

3 

83 

3 

420.41 

258.71 

13.92 

wO M 

3 

84 

3 

278.72 

218.04 

17.67 

226.99 

3 

85 

3 

302.34 

202.48 

26 . 63 

365.00 

3 

86 

3 

377.67 

296.53 

25.80 

149.76 

3 

87 

3 

576.25 

352 .84 

47.54 

224.47 

3 

88 

3 

339.93 

312.38 

123.67 

140.24 

3 

89 

3 

232.04 

198.77 

55.11 

176.29 

3 

90 

3 

464.19 

280.69 

13.79 

157.81 

3 

91 

3 

368.29 

265.54 

26.23 

276.63 

3 

92 

3 

345.92 

244.79 

8.25 

282 . 03 

3 

93 

3 

228 . 37 

210.81 

43.31 

444 . 88 

3 

94 

3 

407.99 

309.23 

61.47 

435.35 

3 

95 

3 

356.67 

221 .71 

11.79 

162.85 

3 

96 

3 

161.12 

120.80 

94.77 

303 . 79 

3 

97 

3 

565.34 

323.23 

84.90 

258.93 

3 

98 

3 

285.99 

250.90 

36.04 

202 . 35 

3 

99 

3 

157.66 

163.92 

135.99 

268 . 44 

3 

100 

3 

199.90 

167.37 

100.30 

123.78 

3 

101 

3 

167.37 

100.30 

7.88 

207 . 1 1 

3 

102 

3 

546.93 

328.59 

71.90 

278.03 

3 

103 

3 

317.65 

224 . 85 

28.33 

137.72 

3 

104 

3 

319.35 

218.13 

43. 12 

108.75 

3 

105 

3 

347.44 

269.90 

31.71 

146.03 

3 

106 

3 

238 . 53 

214. 11 

61.30 

208. 10 

3 

107 

3 

310.32 

184.96 

23.59 

313.26 

3 

108 

3 

256 . 55 

169.08 

8 . 56 

289.64 

3 

109 

3 

178.76 

147.59 

48.88 

526.85 

3 

110 

3 

168.35 

127.57 

42.88 

300.33 

3 

111 

3 

204.34 

151 . 13 

38.07 

509 . 46 

3 

112 

3 

260.70 

219.60 

44.72 

217.50 

3 

113 

3 

254.93 

163.65 

8.60 

252 . 44 

3 

114 

3 

340.71 

224.23 

9.37 

134.95 

3 

115 

3 

357.05 

214. 14 

28.84 

120.29 

3 




•16 

3 

316.46 

500.25 

249 . 77 

357.81 

3 

•17 

3 

336.63 

183.97 

27.22 

339.56 

3 

18 

3 

151.85 

83.20 

302.88 

369 . 36 

2 

19 

3 

334.91 

243.45 

739.22 

126.21 

4 

20 

3 

506.23 

303.24 

18.68 

144.06 

3 

21 

4 

386.41 

257.61 

543.88 

7.51 

4 

22 

4 

358.41 

242.41 

415.54 

8.77 

4 

23 

4 

287.54 

232.03 

442.22 

2.84 

4 

24 

4 

329.14 

253.47 

869.47 

15.25 

4 

25 

4 

319.50 

261 . 4 1 

529.96 

10.36 

4 

26 

4 

317.69 

255.85 

493.23 

0.24 

4 

2.1 

4 

340.51 

286.97 

957.05 

21.82 

4 

28 

4 

279.72 

236.39 

681.81 

19.00 

4 

29 

4 

262.98 

255.69 

362 . 80 

27.30 

4 

30 

4 

260.82 

229.80 

323.59 

22.63 

4 

31 

4 

277.34 

251.29 

597.87 

6. 86 

4 

32 

4 

320.07 

278.70 

470.60 

3.83 

4 

33 

4 

379.81 

269 . 1 1 

744.64 

4.85 

4 

34 

4 

379.20 

264.32 

696.35 

3.96 

4 

35 

4 

252.93 

215.62 

464.23 

12.26 

4 

36 

4 

271.49 

231.45 

176.25 

75.62 

4 

37 

4 

238.17 

195.47 

157.23 

89.96 

4 

38 

4 

312.63 

256.94 

472 . 49 

0.75 

4 

39 

4 

299.32 

265.81 

495.62 

2.07 

4 

40 

4 

337.51 

270.37 

632.69 

3.42 

4 

41 

4 

361.94 

272.68 

683.67 

4.73 

4 

42 

4 

327.07 

239.05 

1147.34 

91.65 

4 

43 

4 

574.72 

266.73 

1115.70 

186.69 

4 

44 

4 

281.62 

262.87 

919.93 

39.41 

4 

45 

4 

329.50 

279.22 

585.49 

7.08 

4 

46 

4 

272.01 

244.05 

504.25 

4.78 

4 

47 

4 

311 .56 

265.12 

638 . 48 

3.45 

4 

48 

4 

305.05 

218.34 

450.36 

9.23 

4 

49 

4 

335.44 

343.43 

531.03 

89.18 

4 

50 

4 

247.57 

243.72 

896.14 

41.01 

4 

51 

4 

420.93 

279.14 

1158.49 

32.32 

4 

52 

4 

327.98 

272.25 

1225.07 

55.10 

4 

53 

4 

387.18 

284 . 42 

423.61 

15.12 

4 

54 

4 

580 . 96 

311.95 

641 .54 

46.18 

4 

55 

4 

474.37 

277.92 

639 . 66 - 

27.92 

4 

56 

4 

265. 16 

231.60 

442.10 

5.71 

4 

57 

4 

253.79 

120. 14 

887.16 

349.51 

2 

58 

4 

352 . 28 

254.56 

538 . 33 

1.45 

4 

59 

4 

323.23 

271.30 

541.16 

5.22 

4 

60 

4 

261.84 

242.93 

664.32 

12.74 

4 
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1 


T-T-T-T-T-T-r-T-T-r-T-T-r-X-T-T-TrT-T^T-t-T-r-T-t-T-T-T-T-T-OJCOfUCUai 


co«t<t*'^cOinin’si'OODjWin 

OO^inC 0 CUNCVl«^V"n 4 -O 


C* 3 > 0 O<JNt 0 ONNin 00 C 0 

cxicuaiDJCxicxicucxjcxjair-cxj 


tocoininNcuo^inc* 3 '^co 

in>o^T-<j£hOT“rov-oco 


O^WDJCUCONCXJDJNOm 

in^oo^w^ocuN’r-tncoin 

OJCUWaiCUCUDJJXltvJCUCUCU 


ocxjN«!?'^ini>T-o<iin 

>OI>‘^CU<JrOC>COCU' 43 l> 


CUCOcOCU’T-cOthONOO 
KtvO'^'=t-^incoNo^^oo 
CUCUDJCUCUCU' r-CU 


intOOONON-sOtO^ChChCONCUinNOiODJOO^ininNOOODI'^O^f-^O^Ch 

Oh«.'^J>N^> 43 N^in^Chin>^NNOJCOCOtO':t<iirtl>COCOOJ>OOOiriin’r ->0 

CO’^Oin<t«tCOCn’ 4 "^^ 0 '^cOcO'O^NT-in' 43000 C 01 >* 4300 ^^COcOCO’ 5 tin 
rOOC>rOChDJOONGC>NNCh<ICVlir!inCh'CO^NCOOOI>ChNO’^“^'^CaNO 
'^CVl'^CU'^CUCUDJ^^ ■r-r~r-T~T-fU‘r-T-T~CU^CUCXlDJ'^'^'^DJCU 


njnT-c^turococxj^cc^c>NO>ocoinNvOc^oroocucooD^incxjin*^^r-oi>N 

Nj>>oo{><i<ichinT-oaDuicO'Oi>NinNc>oinoO' 0 ’^<iincuNoc>«^cucei 

otoroiC>iriN<t‘Ni>cof\iniri!>^(>ja 3 CU'^^cxjcuwi>T-a 30 T-i>N^>ocojrt<i 
■^NChtO^DJCUinOCUNNOinChOCUin-ONT“r-OinCOCh^£>thiri< 3 NCO«:f«D 
in'tin< 5 >GN' 43 NOsOin’^iii>o< 5 iA<>inoinNO^-sO'«o<i*^> 0 ' 4 Din^< 3 '^ T-cuctj 


I 

I OCOCOOtOrOCOCviaiO^<}«=ti»'<3COirOCONcOI>ChNNcOO^*=fCUcOCOCON-43CO 
3 } ^ininNO'NN^‘=t’co^cr)^^ihcuNNOcoo<iococoinNoo£>cO‘^irii>^c *30 
03 1 

I cO’r'inoNinin^cin<iT-T“^inNNcoN^inco(>in<iN<icoNN<icocMrviroco 
{ T-i>cO<tN«CUOOCviC>OT-T-NNO<t>Cr0^tClNCUCUChOC^N<t^CO COOCXl 

I 

I 

( 

! 

{ 

1 

D I C^*:t-^WChCXjcOCOCO>^cOcO'^^NCXJChOI>’=t(hCO*^'=tcO*=t>43cON*^<>‘c-COCiJ<l 

m 1 o^DJin'^‘^c>c^ininoj'Ococ>"^NO^oNNCoco'^’ 5 tcuocu'^i>Ni>":“iriirio 

\ tOCU^NrOCOCO^O'^S>- 43 CxiCUCU^NDitO'^NOL^lcOinN^OJ«^^C>COOON'^ 
j Otai>cONrOCUrOcCn<5NnDJCGl>D"l>iri'^CO^CUCUcOOcOOCOCOD4-^ tON 
I N^tcvit^cucurocxicxJcuc^cunruT-wctjcotxjDiroDJDjnaicucutocxjco 
1 


1 

3 I O^J>’*:“v-i>cOCUChOO£O^CON<lC>inChChOcOCOjri<t>^DJODJC>cOWcO«?tCU 

12 I «:|-CU^O^DJ^U^’^l>rO‘?t‘01’^NN«C'431>ONNC01>*^«5i'(>Ca'^iXS^^<fCOtU 

I CUinCUDiCtjCU*^OCOT-0'^'OCO‘^NtOincONrODiNT“T-|>CONinDJChOcO^ 
ico^ T-T-CUtO*^ T-CXJ^-^COCU -r-Nifix-rs. 

I in N in ^ 


I 

• [ 

a \ 

2 \ 

I 

m ! 

^ I ^^^^^^^^r---r*T-T-'T“”r''r-T-^T-?rT-'r''r-‘C“-r-T--r-T“T~T-TrDJC\jCUDJW 

-J I 
O I 


U 

a 

2L 

m 


f 

I 

I 

1 

I 

I 

I 

I 

\ 


tM to ^ 


in< 5 Ncccho^cuco*=t-inNO 


NCCOO"!-CUC 0 ^in^ 0 NC 01 >O*^DJi? 1 ^iA 

T-T-Tr-cucuwfucucucucuajixifocococototo 
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602.69 4.38 26.07 122.05 15.87 91.09 


7S 


curufUDjcuwcucuDjrowcucxjcxjCuwcufuwcucxjcxiruDJCoc^ncococococorocoixic^c^c^c^c^ 


%0 CO 111 

in N in CO 


N^cochincuo^coincoo 

«;f*r-OCU>0^’r-OC0>0Ch 


CO 0^ £> O CO CO ^ 

N in ^ N CO in 


in'^"tn'€incoO‘^c>'^<ii>cooin'^^ 

inNCUNi>ocO'C^ooNcom<j«t^ 


inOCMN’:“^COt>CU<3CO<tNNONin^£tJ-4DOin^01>'^COCODJ^<l>!}'OO^cOO‘^ChCU 

cONO£>0^'^OcC!NNO^cOcOOcONChcOcONNinCh^^’5t*<lNONONh«.>0<lcOONNCXjNr^ 


Ch^NrNco<iincho^>o 

inNCOCOCUCOcONinr- 


tnT-O^N<3fOl>i>NO 

T“T-ru<t-^co«tco ^ 


cOO^'^cOinOi>NcO 

co<3incooinNOi>T->43 


CO^JNCONOCUOWO 

cxico^cvi^r-cocxicoaiT- 


ooocococoroO"«!to 

OONCr-COCOCOCO**:-*^ 


’T-NOcOONN'ct'-ONin 

T- £0 T~ ^ in 


inNOr-t-COTr^OCO^f^ 

CChONCUcOcOa'O’^ 

co^c^Ncoco^dlnco 
cu cxico^^inm*^ 


0‘?3‘*^NNin'Oi>-43^^r^NNo^^r-incu*^^cooO'<icococooocoin^^tnNO<t‘Nco 

co<ioNinco<tincocoocoN<iCocoN'0*^-oococooDJ-^~orv^co^inch^c>ctj^o*^i> 

^CXJO^OCOO^<lCh*^CUcOCO'^'ONCUOCO<20l^l>in'^'r-T-‘0^cO'^NCOCXjOinOJ*^0‘^"Dj 

aix-'<|T-CONcOCUCU«!t-0^<fO'OCh^NNO^tOOOin>0'^COCOCO'^“NinC>cOcONCOOCOin 

T-x--KrtM-rrrrfXlx- 


cij>ONO^ininchcoinin^ctjinin'^ON’?t‘>ODJCxi*>:“‘^'^ONcococo^coocooco^in*?t-i>^ 

incocoO'r-x-oco<fcoN<tcococuincucocO'«oinoj'^-5-ocufUinin'^coo^NocxjNininco(xi 

T^O^NCONN<t^CUind-odtONNNOCU(>inONOChCOCOCtJCUCODJ^COsOCUl>C^CODJCUCU 
cuT-<t^cu<iinco T-cO’^'^incu’^cocu't^'^cu^coDj co-^^ djcxj 


OT-N‘^cor>.incuin^^N'r-inoincococuDj<tcuoincocoin*^coccoco^T-o<tin<r*!t'0 

coNDiONcochinNO>oin'^r^inco'=tinr^>c3DjNininco«tcoi>in'^NOoinsocxicoDjOinNN 

^DJNN^CO^‘^"0>OcO'=i‘OC^NN’^ini>Oi><3NCU'^N«5tCOOcODJNN'OCxJCONCOC>CU 
T-CO CU‘r“O-«r“T-tUCUC0-^’^*^CUI>'r-T“C0’^ 


^<5incOr^^<t‘OJN*^cOcGCON'OincONCOa'COCUinCUN<jiCO^N^CO 

cocuNNincu-^^rO'43N'Ci>c>in'Ococuin*^!>cu*?“^^NCu<3ChCuco 


<r o ^ to 

o m o *= 1 - ^ <3 


't 

N CO 


NO^O^ODJ-OC^C^NOCO^^OJ<3^0^COlncOCOcOCOC^JC^■^^'^COCOC^CUC^'-OC^>^tlnC^C^T“ 

china3'^~<i*'”'^^chcoco^coin’^-ocococuco^<3co<onin*^incucootocoo^(>^co<i^ 

inNO’?t^in*^^inch<J'5tin<3<’N<3^inN-^>oi>>ocuccNoi>0’^Ni>cofvf>v.i>T-inin«o 


CUCUCUOJCUCUCtJCUCtjCUOJCVlCUaiCUCUCUCUCXjCgOjDJCUCUCUCOCOCOCOCOCOCOCOCOCOCOCOCOCOCOCO 


'ONcOO^O'^CXiCO 

cococoro*^*=t<“*'t 


^in<3NCOD^OT-CXJtO'=^tn-^NCOO^O*^CUCO^in<INcO!>OT*DJCO«^in-43 

^ ^ ^ ^ ^ ^ 0^ tf^i 11^ trt tj^ wTl .jn «J*i ^ vH ^ N.. N.. JV, 



492.89 5.93 32.39 140.05 36.52 96.02 
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OOlXiCXicOCOCO’^Cij'<!<JOONiritriCUNOcG^iri^cO^>C'^NcOOCUCOCh'^’^~O^OJNN'^" 

cO'^CrJcOCO<5>CC>>€il1-?“NCOWC^WO^Ch^DJ'ttl>WcODJcO"t<lin'?t'=t'<linrO|vcONinN!> 

cur^oinDJcocuocuococ^'COin*o^o^cho^O‘^o^co'!^cu^inin^N"^iriitQN^^T--sODitri 

ixiin^<iNO>oirj-sOiri^inNocxj’r-rN.coN«5f^cuinD3conNDjN<ic*7cocGND3Noo<jcoin 

T- -r- T- T- ^ r~ T- r- r- T“ T- ^ rr -r ^ tli r- CU 


o^Noc^c^co<l^oco«^coo^NOlrlcc)lnu1’^tO'ONlnl’5f^^<;^<tDJcoN<JaJ‘^“OtO'43c^lncoco 

o-<iJ>coococo':tcooi>i>oNCh>0"^<JNcoochWcocoN^iricuNrocu-^*^*^’^inNOino 
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Table 5.6 Results for MSS Test Data 


Number of samples correct-l y classsified as 


Predicted Water 
class 


Urban Forest/ Uncultivated Accuracy 
Veg. X 


Water 

Urban 

Forest/Veg . 
Uncultivated 


1007. 

1007. 

957 

97.5% 


Table 5.7 Results for TM Test Data 


Number of Samples Cor roctly Classified as 


Predicted Water 
class 


Urban Sub- Forest Veg. Unculti- Accuracy 


urban 


vated 


Water 

Urban 

Sub— urban 

Forest 

Vegetation 

Uncultivated 


— 

— 

— 

— 

— 

1007. 

28 

02 

— 

— 

— 

93.3% 

02 

28 

— 

— 

— 

93.3% 

...» 

— 

27 

03 

— 

90% 


02 

— 

28 

— 

93.3% 


- 


01 

29 

96.6% 
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As can be expected 3 samples from the predicted forest 
class are classified as vegetation due to lack of dense cover. 
Two samples from vegetation classified into sub— urban class 
suggest the increase in vegetated areas. 

/ 

For a classifier increase in the number of classes is 
associated with decrease in the classification accuracy. But 
the analysis using TM data indicates that even after 
increasing the classes from four to six, the classification 
accuracy is maintained over 90 percent for all the classes. 
The Survey of India toposheet used for taking training set 
data has 1:50,000 scale. This has resulted in taking the 
training data accurately and a correct cross check with the 
test sample data. 


Table (5.3) for the training data set for TM shows the 
values of tranformed vegetation index against each class. It 
can be seen that the value of this vegetation index increases 
with increase in vegetative component. For the present 
training set a value of more than 70 has good vegetation 
content. For forest, the values are still higher as shown in 
the Table. This vegetation index can be used to reclassify 
the vegetation types. In an area which is intensely cropped 
and has various types , it is possible to classify it in 


different crop zones 
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Although six categories have been taken in TM, the 
features such as roads, canal etc. fall in one or other class. 
Such features can be called as micro clases. Roads, if 
metalled would come under the urban and sub— urban 
classes. Typical village roads and earthen roads would 
be classified under the bare soil due to close resemblence 
in reflectance. The use of classiifier is limited in such 
situations and use of imageries should be preferred. 


5. 4, Resulbs of* Image Processing — 

The simplest way to get the image display is through the 
computer line printer. One band density slicing is done to 
slice the entire range of reflectance values in given number 
of classes in the scene. Fig. (5.1) shows the line printer 
map for band-4 data. It can be seen that except the river 
pattern, it does not convey substantial information. It 
essentially lacks the quality to create an impression on the 
brain. But the method is very simple and easy to apply since 
no special system is needed. The display is very quick. 

B.S Coii^arlson of Map and the Imagery 

To prepare a map is a hard and strenuous Job. It 
requires labour and heavy expenses. It is difficult to take up 
the exercise of updating the map at very short intervals. An 
imagery on the other hand is a data product generated from the 
digitized data. Although the imagery does not give the 
information in complete details^ it can be used to study the 
temporal changes. 
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Figures (5.2) and (5.3) show the details covered in the 
map and the imagery respectively -for the Lucknow area. The 
two figures resemble quite closely in respect of the main 
features. Particularly the river pattern of Gomti compares 
very well. This suggests that a time lag of ten years has not 
changed the river course. 

Figures (5.4) and (5.5) show the comparison of the map 
details and the details covered by imagery for the pattern and 


flood plain 

of 

Ganga. 

. The general 

pattern 

agrees in the 

two 

figures, but 

within 

the natural 

limits 

the 

river has 

a 

tendency to 


shift 

as shown by 

the imagery 

. There 

is 

continuous 

change due to silting 

and 

deposition 

of 


sedimentsr typical of anastomatic pattern. 

5.6 Con^arlson of MSS and TM Digital Images - 

Digital images are very useful where the facilities 
exist to carry out image processing techniques. Apart from 
the simple one band digital images, filtered and zoomed images 
are used to enhance and enlarge specific features. 

Plate 1 and Plate H show the digital images in band 2 of 
MSS and TM. Forest and vegetative features are highlighted. 
The area covered by the TM image is about one seventh of the 
coverage in MSS. It can be seen that for studying the 
lineaments MSS image is more suitable. 
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Fig. 5.5 Details of Imagery (Flood Plain of Ganga) 



I 
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Plate 1 



Plate 2 
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Plate 3 and Plate 4 show the -false colour images in band 
2. False colour images are useful for differentiation when 
difference of grey level is difficult Itto perceive. 

Plate 5 and Plate 6 show the digital images in band 4. 
In KSS image there is little variation in level in surrounding 
area around urban land. There is greater coverage in NSS and 
effect curvature of the screen is also pronounced. The TM 
image shows two bridges over the river Gomti. 

Plate 7 is the false coieiU' image in band 4. Urban area 
near Charbagh is misclassif ied. Brown colour shows the sub- 
urban fringe. It also covers the land recently opened up for 
built-up activity. 
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Plate 4 
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Plate 5 



Plate 6 
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Plate 7 


CHAPTER VI 


CONCLUSIONS AND FUTURE RECOMMENDATIONS 
6.1 Conclusions 

The present study establishes the utility of various data 
products for the analysis of the remotely sensed data. The 
two modes of data collection i-e. mult ispect ral scanner and 
thematic mapper have their limitations. 

By supervised classification for the samples we arrive at 
the following conclusions: 

(i) For general analysis MSS data can be used but it can not 
classify the micro classes like roads, railway lines etc. due 
to coarse resolution of 79m. 

(ii) TM data is very useful when the classes in the scene are 
many. It is possible to classify the microclasses. 

<iii> Handling of TM data is little difficult because of the 
volume of data. Also it costs about eight times the c ost o f 
the MBS data. Hence for general purposes MSS data shall find 
its use. 

The use of imagery alone is not very useful unless 
available with the digital data. This is due to the fact that 
the digital data can be controlled in a desired way but 
imagery can not. It is very easy to produce the digital 
image for selected area and perform other techniques of image 
processsing . 
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6. S Future Recommendatl ons 

Bayes classifier is very efficient classifier. In an 
area where there is wealth of variety, it is a very useful 
tool when TM data is used. To get the best results from the 
classifier all the seven bands can be used. Different 4 band 
combination can also be tried to get the best combination of 
bands. This would save the computational time. 

For complete classif iat ion of the image data vegetation 
can be re-classified into sub-classes using vegetation 
indices. For the detailed study of the urban area, image data 
can be zoomed and then contrast stretched. This would 
differentiate amongst residential, commercial and industrial 


zones 
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PROGRAM DENSTY.FOR 


THIS IS AN' IMTERACTI'.'E PROGRAM FOR MAPPING THE 
SURFACE FEATUURES BASED ON THE DENSITY SLICING METHOD 
PROGRAM IS FOR EFFECTI','E FOR 10 CLASSES 
THERE IS NO LIMIT ON THE LENGTH OF THE INPUT BUT 
THE NUMBER OF PIXELS PER LINE SHOULD NOT EXCEED 96Q 


INTEGER IPIX ( > r INTMAT {96Q) t GREY ( 10 ) 

INTEGER CONTOR(10)7lNT'.'L<10> jUPLIMdO) rLOWLIMdO) 
0PEN<UNIT=257DE'.'ICE=--'DSK' rFILE-^dNPUTd 
OPEN < UNIT~£ y DE'JICE«' DSK ' y FILE- ' MAPI ' > 


0PEN(UNIT=:7yDE'^ICE~'DSKdFI'^E .-'MAP2' > 
OPEN ( UNIT-Sy DE'v'ICE- ' DSK ' y FILE*- ' MAPS ' ) 
OPEN ( UNITES y DE'^ICE=:" DSK ' y FILE-" MAP^ " > 
0PEN<UNIT=*10yDE','ICE--'DSK' yFILE«"MAP5" ) 
GPEN<'JNIT-llyDE',»ICE«'DSK' yFILE™"MAP£' ) 
0PEN<UNIT«-42yDE','ICE-'DSK" yFILE=^"MAP7' )' 
0PEN<UNIT-13yDE'^ICE-."DSK" yFILE-"MAP8' ) 


C 

c 

10 


20 


C 

C 

C 


^0 


•jr> 

w w 


cr,r, 

W vr 


INTERACTION BEGINS 


TVOC I A 

i < • *» AW 

FORMAT(" ALL TERMINAL INPUT IS FORMAT FREE"y//y 
ISXy'TYPE IN THE NUMBER OF CLASSES FOR SLICING' > 
ACCEPT -YyNCLASS 
TYPE 20 


FORMAT<SXy 'TYPE IN CHARACTERS FOR REPRESENTATION 
1 OF THE CLASSES') 

ACCEPT 25y (GREYd) yI=lyNCLASS> 

FORMAT dOAl) 


TO GET THE CONTOUR LE','EL AND INTER','AL OF EACH CLAES 
FROM TERMINAL 


DO 30 I-lyNCLASS 
TYPE "10 y I 

FORMAT'SXy 'TYPE IN THE CONTOUR LE'JEL AND INTER'.'AL' y 12) 
ACCEPT .t. y CONTOR ( I ) r INT'JL < I ) 

I ic'i TM / T \ ^r«n>.iTnn< / T \ xTxiTiK 

LOlv'LIM < I )-CONTOR ( I )-INT'.'L < I ) 

CONTINUE 
TYPE 50 

FORMATiSXy'TYPE IN THE LENGTH AND BREADTH OF INPUT 
ACCEPT Y.yNLINEyNPIX 
TYF’E 55 

FORMAT (SXy 'TYPE IN THE NUMBER OF OUPUT FILES') 

ACCEPT -t-yNOUTFL 


IIT % 


c 

C INTERACTION ENDS* 

C T£) INITIALISE THE INTERMEDIATE MATRIX INTMAT 

C 

DO 65 I=lyNPIX 



65 

C 

C 

r 


6Q 


70 

C 

c 

c 


TG REiM! AND CLASSIFY THE IHPUT DATA 
ILOOP^O 

rnKiT Ti<u ip 

UrWl%i iABtwiL. 

FbEAB {25ySi > (IPIX<I> 7l"“l rMPIX) 

Ti nnD«-.T| nno jl i 

^iwwwi I «k 

DO 70 1**^1 jNF* IX 
DO 70 J=irHCLASS 

IF< (IPIX<I> .LE,UPLIH( J>) .AHr,<IPIX<I> .GE,L0WLIM<J>> 
1IK'TXAT<I)«GREY<J) 

CONTINUE 

Tc trcnsfcr the rcculte frcni INT.MAT to outrut filcc 

•f” r 


icr the line rrintcr 


ulwi 


OA 


nn OA Tssi -KiniiTci 

A.* imt WB/ Im 

ypi TE < NUNIT 7 90 > < IHTMAT < J > 7 J”M - N > 
FORWAT<120A1> 

Ml IKITTc-.MI IMTTJL 1 

M 1 1 

N«MF119 

CONTINUE 

IF (ILOOP,LT,NLIN£) GOTO 60 
STOP 

CKin 
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READ. FOR 

PROGRAMME TO REAAD CCT ON ND-560 COMPUTER 
PROGRAMMED BY K.V RAO 
INTEGERifI IBUF(4000) 

OPEN ( UNIT=ED , FILE= ’ INPUT : DAT * ) 

DO FOR K=1,2 

ISTAT=MAGTP ( 10B , IDUM , 40B , IDUM , IDUM ) 

END DO 

DO FOR 1=1,1680 

ISTAT=MAGTP( 16B,IDUM,40B, IDUM, IDUM) 

END DO 
DO 3 M=1,290 

I ST AT=MA6TP ( 26B , IBUF , 40B , 4000 , IDUM ) 
WRITE(25,4) (IBUF<J) , J=1300, 1589) 

4 FORMAT (5014) 

DO FOR N=1,6 

1ST AT=MAGTP ( 16B , IDUM , 40B , IDUM , IDUM ) 

END DO 

3 CONTINUE 

CL0SE<25) 

STOP 

END 



C prc:2i'oni RECORD. FOR 

C 

C 


r 


C 

C 

c 


PROGRAM TO/:EAD AM IMAGE LIME OF THE CCT 
THIS PROGRAM READS OME RECORD OF CCT, THIS IS VALID 
IMSERTIOM OF 5 r lODUHMY BLANKS OR AFTER CONVERSION 
CCT FORMAT INTO DIC-10 COMPATIBLE, 


AFTER 

OF 


C 






* 4 ^ 4 ^ 4 ^ 4 ^ 4 ^ 4 ^ 4 ^ 4 ^ 4 * 4 ^ 4 ^ 4 ^ 4 ^ 4 ^ 4 ^ 4 ^ 4 4 ^ 4 ^ 4 ^ 4 ^ 4 ^ 4 ^ 


11 

1200 


21 

101 

10 


INTEGER OCT ( 1 125 ) 7 BTl r ET2 7 BT2 7 BT^ r BT5 
OPEN < UNIT--53 7 DEVICE=.' DSK,' > 

OPEN ( UNIT-50 7 DEVICES' DSK ' ) 

OPEN < UN I T--! 5 - DEV ICE-i' DSK ' > 

0PEN(UNIT--20tDEVICE-'MTA1 ' 7 mode- 'DUMP' 7 RECORD SIZE-11257DENSITY 
1 -'1£00') 


N0--3 

READ (20) OCT 
DO 1200 1=1 7 1125 
N0=N0+>1 

BT1=0CT<I>/2S:^2S 


BT2= < OCT ( I ) -BT 1 J^O-t-tOS ) /2S;S'.20 

BT3= ( OCT( I >-BTl.Y2!{'..t.29-BT2:}:2S:S’.20 )/2*S:12 

CCT ( I ) '“BT1^2*T*#2S“'ET2A2M-S^20“-^T3*^2.'i'*X 12) /2if.# ^ 
WRI TE ( 53 7 1 1 ) NO 7 BTl 7 BT2 7 BTsTbT^ 

FORMAT <51 5) 

CONTINUE 


CALL CONE 

CL0SE(UNIT=--207DEVICE='MTA1' 7M0DE.-'DUMP' 7 RECORD EIZE--II 257 DEN 
1SITY=:'1600' ) 

STOP 


END 

SUBROUTINE CONE 
DIMENSION NAT<25> 

COMMON NAT 

OCI.it MTi ^.1 

* W w 

DO 101 I=l7lS0 
READ < 53 7 1 0 ) ( NAT < K > 7 K-1 7 5 ) 
READ < 53 7 1 0 > < NAT ( K > 7 K=6 7 1 0 ) 
READ ( 53 7 1 0 > < NAT ( K ) 7 K-=l 1 7 1 5 ) 
READ < 53 7 1 0 > < NAT < K > 7 K=16 7 20 > 
READ < 53 7 1 0 ) ( NAT < K ) 7 K=21 7 25 > 
WRITE (50 7 21) NAT 
FORMAT <25 15) 

CONTINUE 
FORMAT <5 15) 


F^ETURN 

END 



n n n n n ti o rs o 


program: 


PIXEL. FOR 


Thic r-rorircni iz to rood o record of the CCT to det the 
drou level of cortein rijsel whio hoc to be fed into 
rredrom uhilc under cjjecution. 

This rrodroni oloc rrovideo the reflectonce voluoo 
of five rixelo on either side of the wonted one. 


30 


i n/\r\ 

a. w w 


31 


ill 

U. JL. .1 

i 1 2 


I NTEGER CCT < 1 1 25 > i ET 1 7 BT2 ? BT3 r BT >! » E YTE < 1 000 > r ABC < 1 1 ) 

INTEGER RECORD 7 BAND 

OPEN < UN I T»=^53 r DE'.' I CE** ' DSK ^ ) 

0PEN(UNIT=507DE','ICE«'DSK" xFILE-'PIXEL' ) 

OPEN ( UNIT=20 7 DEVICE^' MTAl " 7 MODE- ' DUMP ' r RECORD SIZE=1 125 7 
1DENSITY='1£00') 

TYPE 30 

F0RMAT(" TYPE IN THE PIXEL NUMBER PLEASE') 

ACCEPT S:7lPIX 
READ < 20) OCT 


J=1 

DO 1200 1=^1 7 1125 

OTi =.nf'T r T ^ /'i'if'lf'iO 

BYTE<J)=BT1 7 J=J+1 

BT2- ( OCT ( I ) -ET 1 ;i:2.t.S:2S ) /2-'*:^i20 


BYTECJ>-ET27 J-J+1 

ET3^( CCT ( I ) “ET1S'*2^-^29“~ET2m'-2XS'-20) /2S'-S^12 
EYTE(J)^ET37 J~JT1 


p ^ OCT ( I ) — BTIS 2 Xm^ 26— ET2S^2M^-i^20*“BT3i2S*S*12 ) 

BYTE< J>=BT'l? J=J+1 
K- 1 * 1 7 I P I XU- 1 P I XT ^.00 


IF (K.GE.IPIX'J) GOTO 31 
CONTINUE 

CLOSE ( UNIT~20 7 DE'.'ICE- ' MTAl ' 7 M0DE= 'DUMP ' r RECORD SIZE'-l 125 7 
1DENSITY--='1600' ) 

RECORD-BYTE < 7 ) -f BYTE < 9 ) ^ 2'56 
IBUFF*f:ECORD/'i 
P.AND=RECOF:D-IBUFFS:^ 

LINE'=EYTE( 11 )+BYTE(12)J}:25£ 

WRITE<507lll) LINE 7 IPIX 
F0RMAT<15X7 'LINE NUBER- ' 71 ^. 7 ' PIXEL 
WR I TE < 50 7 1 12 ) BAND 7 RECORD 
FORMAT ( / 7 15X 7 ' BAND •■=-' 71 ^ 7 ' RECORD 
I ZERG"*^ 1 

TO COUNT THE LENGTH OF THE INITIAL ZERO FILL 


NUMBER -/7l5) 
NUMBER « '715) 


i O*^ 


^ W 


113 


1*0 125 I^lZr‘100 

IF <EYTE< I> .EG.O) IZERG---IZER0+1 
TYPE i'.7lZER0 7lZER0 
IPIXEL’^-IPIX+IZERO 
IPIXU-IPIXEL + 5 
IPIXL*=IPIXEL - 5 
J=*l 

DO 113 I^IPIXLrlPIXU 
ABC<J)=BYTE(I> 

CONTINUE 



ll'i 

115 


WRITE (50 r 11 1 } AEC 

FORMAT <// 7 15X7 'REFLECTANCE COUNTS ARE'/Zt ISX r 11 < 15 ) ) 

UDTTr -1 1 «=:\ 

F0RMAT</7>1'1X7"''7/x1'1X7 ' ! ' 7 / 736 X 7 
1 ' WANTED PIXEL ' 7 // 7 15X 7 55 ( 1 H- ) ) 

TYPE ♦xLINErlPIXrBYTEdPIXEL) 

STOP 

END 



c 

r 

r 


■* '^00 


C 


THIS PROGRAM READS A RECORD OF THE CCT 
OUTPUT UILL BE ONLY A PART OF IT. 


BUT 


THE 


I HTEGER OCT < 1 125 > r BT 1 y ET2 r ET2 y BT'J y BYTE ( 1000 > y LOWER y UPER 
OPEN ( UNIT=53 y DEVICE* ' BSK ' > 

GPEN(UNIT=50yDEVICE='DSK' yFILE*'PART' > 

OPEN ( UNIT=20 y DE*,'ICE=' MTAl ' y MODE*' DUMP ' y RECORD SI2E-1 1 25 r 
iIiENSITY='1600' ) 

READ(20)0CT . .. ..... ... 

L0UER=1701 

UPER“-1S00 

J*1 

DO 1200 1*1 y 1125 

BTl*‘0CT<I>/2.t.J}:2S 

EYTE<J)«=BTly J*J+1 

BT2*( OCT < I ) -BT lS:2t.>f'.2S ) /2^.*20 

BYTE( J>*ET2J J-J+1 

BT2* <OCT < I > -BTl.Y;2.’f:-t2S-BT2S:2.l:.'i:20 )/2-t.t.l2 
BYTE(J)--ET2y J-J+1 

BT1-- <OCT < I ) -BTl-T•2^}:S'•29-BT2S’•2^}:^•20-ET3t.2.t.tl2 > /2t-t 1 

BYTE<J>*BT1y J»J+1 

CONTIMUE 

CLOSE < UN IT*20y DEVICE* 'MTAl ' y MODE* 'DUMP' y RECORD SIZE=1 125 
1 DENSITY*' 1600' > 

IZER0*0 

TO WRITE THE DATA 
DO 3000 I^lZylOO 

IF (BYTE(I) .ED.O) I2ER0*IZER0 + 1 
LOUER=LOUER + IZERO 
UPER=UPER + IZERO 

URITE<50tm > (BYTE< K> yK*LGUERyUPERy5) 

STOP 

END 


3000 



20 


C 

c 


100 

r 

C 


C 


200 


riMENSION' X(200) rY(200> rDISCSO) jCLSTXC( 50> 7CLSTYC<50 
INTEGER T0GC200> rNCLST<50> 

IiATO. NCLUSTt NPTS 7THRESH/'! 7 2007 1,25/ 

OPEN< UNIT-31 7 PILE-' N, ' ) 

WRITE ( 22 7 1 1 ) NCLUST 7 MPTS 7 THRESH 


CONTINUE 
DO S'! I=:l7lO 


NCL£T(I>-0 

OpUTKin Ti 

Wat. 

II =1 
NCL =1 

READ THE FIRST SAMPLE VALUE 


READ(3l7S:)X<II)7Y(II> 

ri CTVr/MPI \ — V/TT% 

m CTvr/Kiri % 

WWW* tw%i%ww^ '^i 

TAG < II) =NCL 


NCLSKNCL) -NCLST(NCL) + 1 


READ THE NEXT SAMPLE VALUE 
CONTINUE 
11=11+1 


READ(3l7.t>X<II)7Y<II> 

CALCULATE THE SGARE 
DO 100 1=1 7 NCL 

DIS < I ) =DIST < CLET’^C ( I ' - oxyr / t\ .v 7 tt\.y/tt>% 

pnXITT Ml ip 

Wrwi* « Atf.wiw 


FIND THE NEAREST CLUSTER CENTRE 


riTCMTIil — riTC 7 1 % 

JU«.t.Wtt.*.l« '“**.♦ a. W % .k X 

Tip OAA T-i «Mri 

A.*W WWW .*.^.*.^t*WW 

Tr/riTc/T\ pr riTowTxn pn Tn nAp 

4*>i \a.*xw%u. / irwi vw.t.witA**# WW (W WWW 

riTCMTM — riTC r T > 

NCUR 

pnMTTMl IC* 

Wwit • .kt^ww 


310 

c 


TEST UHETHER SAMPLE 

MTAC-rcT n I icTcp 

4 \W««**t— W ( WWWW I wr% 


IS WITHIN THRESHOLD OF THE 


IF <DI£MIN,LT,THRESH)G0 TO 310 

m CTvr/Kiri \ — v7tt% 

wwwl/«w\fvv/w^ f \ \ ^ X f 

CLSTYC(MCL) 

Mr* —Mm XI 

*\WW — ttww I A 

TAr;/TT% —Mm 

a t iw « .k 4' *^1 a.ww 

Kim er/Kim Mm ct/mpi 

a*wwwt *a\ww^ •— 1 %www* kiawwa au. 

Tp/Kiri nr kipi iieT% nn th iia 

J.I aiawwvwtartkwwwwi^ ww »W www 

po Tfi *=:aa 

W w t w w w 

TAP / T T ^ — KIPI ID 

ittW\.kA# t * WWf * 

AMPf C? —KIPI CT^KIPI1D\ 

i«4%WL>W *~l«Wwwt \««ww»%# 

UPDATE THE VALUE OF THE CLUSTER CENTRE 



r;AA 


C 

c 

330 

£00 


■7AA 

/ v» vr 


-150 

•30A 


^00 

61 

•j 1 

w ^ 

^.1 

>=:i 

w ^ 
1 1 


on 


01 


r*. ^ n 

CLSTXC<HCL'R> =<ANCLS5f:CLETXC( NCUR) iX(II> )/(AK'CLS+l ,0) 

CLST YC ( K'CUR > = ( ftNCLESiCLeiYC < NCUR > +Y < 1 1 ) ) / < <^MCLE-f 1 , 0 > 

NCLET(NCUR) ^MCLET(NCUR) +1 

CONTINUE 

IF(II,LT.NPT£) GO TO 10 
GO TO £00 

IF THE NUMBER OF THE CLUSTERS 
THRESH-=THREEHY.l .2 
GO TO 20 

nnMTTKU icr 

WW«\« 

t.iDTTcr / nn ^ a-\ \ 
nn “7AA T — i « OAA 

WR I TE < 22 ; 51 > 1 7 X < I > » Y ( I ) ; CLETXC ( TAG (ID? CLSTYC ( TAG ( I > > 7 TAG ( I ) 

CONTINUE 

URITE<22,-92) 

nn A «=:a i szt - MAi 

WRITE <227 91 >L»CLSTXC<L)7CLSTYC(L)rNCLET(L> 

CONTINUE 

WRITE < 22 7 31 > NCL 7 THRESH 

DO ^00 I=l7NCL 

WRITE <227£1>I7NCLST<I> 

pOKITTMl ter 

CALL F'LCT ( NPTS 7 NCL? X 7 Y 7 TAG > 

FORMAT (/7 10 X 7 'NUMBER OF POINTS IN CLUSTER ' 7 13 7 ' IS “ ' 7 I'!) 
FORMAT(///rlOX7 'NUMBER OF CLUSTER INTHE BATA AFTER MODE 
1 SEEKING = ' 7l27//7lOX7'THE THRESHOLD DISTANCE - '7F12,£) 
FGRMAT<//7SX7 70(lH-)7//7l0X7'S3nirlc No»'73X7' X ' 73 X 7 

I'Y ' 7 SX 7 'CLUSTER CENTRE ' 73 X 7 'CLUSTER' 7 // 7 SX 770 < 1H-) ) 

CrnDMAT / i OV - T -«s:v , CO O.-IV.CO 0-7V.C0 -I.TV^CQ ■7,c:v.T/>N 

t Wt%««iai 4 / 

F0RfiAT(///7l5X7 'RESULTS MODE SEEKING ^.NALYSIS' 7/7SX770( IH. 
1 7 lOX 7 ' NUMBER 

1 OF CLUSTERS ALLOUED - ' 71 ^. 7 / 710 X 7 / 710 X 7 ' TOT.-^.L NUMBERS OF 
1 SAMPLES = ' 71 ^. 7 / 717 X 7 'THRESHOLD DISTANCE ' 7 FI 2 .S 7 /) 

enpw AT / / / « OV *. -7 A r 1 U- % « / ' rCKITpC / . / - 1 AV • ^ P I I ICTCD ' • TV ^ ' 

• * 0 * f #ww#^# wdi.liCt 0 0 T JL %0 0% 0 wk..wwlk«.i\ #/«i# 

1V V / « Z V . ^ MpMPCrpCUTD/ « /% 

Jk 0% 0 x^0% 0 . /0/«# tiC.j.IT'.ii.AwnAs 0 0 0 

CnCMAT r 1 AV • T , C1 n Z-^V-C10 Z-OV.T.^^ 

{ \ 0k \t 0\ 0 ^ t 0 0 % Jk 0tm. ^ \J 0 » 0% 0 * M^^0%j0Su0t0X 9 0 

STOP 

CALL UERTST 

nAI I IICWMMV 

W « I kwiii. WW» f/> 

Oviri 


c 

r 

c 

r 


r 

w 

c 

r 

w 

c 

C 


REAL FUNCTION DISKXl r Y1 7X27Y2> 

TUTc ciiwrTTnM rAi pin atcc tlic ciiri tticam riTCTAKiri: 

TDCTUCCM Tl.in OOTMT^^ 

A.**— » .<rrW . wA..»iw 

SIk' VI/ V** Vt^ VL' VI* VI* VP VP VP VI* VP VP Vti* VP VP VP V^ VP VP VP VP V/ VP V^ VP VP VP VP VP VP VP VP VP VP VP VP ^P VP VP vP 

riTOT — enc'T / / VI — VO V ^vpo x / y* — vo ^ ’i'VfO \ 

A.‘ ^ w I w t \ I V % a \ Z % jh« a *- . % 9 0k 9 ^ 0 0 

RETURN 

CKin 

^9 % Am* 

VPVPVPVPVPU'^VPV/<*^VPVP^VP>^VP^VP^VP^^VP^VPVP>PVPVP^iIP4'VPVP^VP^VPVP'4’’VVV^]^]^^^^^}^ 

SUBROUT I NE PLOT C NPTS 7 NCL 7 X 7 Y r TAG ) 

THIS SUBROUTINE PLOTS THE POINTS IN A TWO DIMENSIONAL PLANE 



c 

c 

c 


77 


100 


■7rtr> 

/ SJt 


WITH THE APPROPRIATE 
USPLX - IS A LIBRARY 


TAGS DENOTING THE CLUSTER MEMBESHIP 
ROUTINE USED FOR PLOTTING A GRAPH 


D TMCTMCTnM V^MC•TC^ 
TKTTirpro TAr:/MOTc\ 


• V /XfDTC N 

/ 1 «t«* «w/ 

. A / 1 AA\ 

• t * 


VV/OAA. A \ ^ TMAr'4 / •=*. 1 *=> 1 ^ 


FORMAT < 1 'I 'lAl ) 

Tin 1 A A 1—1 -wri 

JLww w 

DO 100 I-lyNPTS 
YYdrJ) -0.0 
DO 700 J -ItNCL 
DO 700 I -IrNPTS 

TC/TAR/'TS cn !\ yw / T . I\ — Y f T \ 


CONTINUE 


I A -NPTS 
INC “1 


CALL U£PLH(X; 
TYPE t7 lER 
RETURN 

CMri 


vv _ Mc*Tc . /» . T yip . 

• f <*%i iw# i>a.atw^ 


IArA;IMAG''.7 


lER) 



A A 


50 


■7 A 1 


TTAO 


•7A:2 


305 


pnXIT T KII 1 C 

WW*^* 

BD 50 I-l7^ 

nn =:a t-t . 1 

ww w 

UVAR < I ) -U',‘AR ( I > ■{- < < W ( J 7 1 ) -y MEAN ( I ) ) .t-tE ) /N 
B VAR < I > -E VAR < I > + ( < B ( J 7 I > -BMEAN ( I ) > > /N 
WAR ( I > - WAR < I > + < < V ( J 7 1 ) -'.'.MEAN ( I > > ) /M 
RVAR ( I ) -RVAR < I ) + ( < R ( J 7 1 ) -RMEAN ( I > > tt2 > /N 

roMTTKn It: 

WW«« • 

i.ic»TTC / on. •T'A-i \ 

• lb. wVa f 

FORM.AT ( lOX 7 ' FEATURE t U'ATER ' ) 

l.ip T TCr 7 00 . -7-7-7 \ 

crncMAT ^ / % 

MOTTC / on . ”2 AO \ 

irti <k tU • 4A.4L. / ww«_.# 

CrnpMAT/ 1 AV.TZ / / \ \ 

UCTTC /OO ^ -7 AX \ 

FORMAT '< 1 2 X 7 ' BAMD'" 7 5X7 ' MEAN ' 7 lYx 7 ' VA.R lANCE ' 7 2 X 

W DTTC /OO . TOO \ 

I » «•> j 4 . \ *. Jw # w V «. # 

EiO 30*1 1*^-1 7 ^ 

WRITE ( 22 7 305 > 1 7 WMEAN ( I > 7 W'.'AR ( I > 

FORMAT ( 1 ^X 7 12 7 2 < 3X 7 F 1 0 , 'J > > 
rnbiTTMi ir 

wwi^* <i*.»4W(b. 

up TTC / 00 .XAO \ 

49 t % * H-m % 4m- 4km. 4 %.r %4 JLm 4 


306 


F0RMAT< 10 X 7 'FEATURE: BUILT UP AREA') 


•7(0-7 


•7AO 


309 


i.iOTTtr 7 00 . •7no\ 

<« 1 >« ^ < 4 . « 4 U. \ 4U. 4mm 4 WW*.# 

yOT TC / 00 . XAX \ 

% X 4 W 4 . 7 W# 

UFbITE < 22 7 202 ) 
nn XA *7 r,-^im.A 

i.iOTTC i nn - XA«=; n t . nMCAKi/- t \ , ©i iac* / T ^ 

WRITE<227777) 

(-‘nMTTMi icr 

wwt> i 

URITE<227302) 

WRITE ( 22 7 30S ) 

FORMAT C 10 X 7 'FEATURE: ','EGETATION' ) 

UDTTP 7 00 .- 7 - 7-7 \ 

MO TTC / 00 .XAO V 

i9t % 4L « %4m,*m.4%mr%44Lm.4 

upTTc / nn . XAX \ 

irrA^*. i fc- % 4m 4m 4 %m4 \4 %mt 4 

Up TTC / 00 . XAO % 

DO 309 I~‘i 7 *1 

WRITE < 22 7 305 > I 7 '.'MEAN (1)7 V'.'AR ( I ) 
CONTINUE 


I.IC'T TCr 7 00 . TAO \ 

99 1 % 4k 9 tm. % 4m 4mm 4 %.4 %4 ^ 4 

up T TC r nn . x i A % 

99 k % 4. 9 9mm \ 4Lm 4mm 4 %4 -4. %4 4 

XI A cnOMAT" / *• Av . ' ccATi IOC ♦ onpK*- / % 

UD T TC / OO - XXX \ 
ypT TC / OO - XAO \ 

99l%m9. 9 9mm %4mm4m.4'9m4%44m4 

W OTT^C ^oo. xax% 

«<I.X a (b. \ 4mm, 4m 4 9m4 %4 %m4 4 

UD T TC / nn m XAO \ 

99 9% ,4. « kmm % 4m 4m / 9m4 %4 Jkm 4 

nn xii 

W V -7 Jb dfc 4 b 4 b # » 



o o o o n n o 


WRITE ( 22 r 305 > I r RMEO.W CD? RVO.R < I ) 

rOKIT TKU IC 


CALCULATION OF COVARIANCES BETWEEN VARIOUS WAVE BANDS 
OF VARIOUS FEATURES 

nn XA 1—1 .iA .... I . I 

cowl W2.T.C0W1U2+ ( < W < J J 1 > -U'NEAN <l))S:(W<Jr2) -WNEAN < 2 ) ) ) /N 
COW 1 U3-C0W 1 W3+ ( < W ( J y 1 ) -U'MEAN < 1 ) > Me < U < J r 3 ) - WNEAN < 3 ) ) ) /N 
CGWl Wl-COWl W^.+ ( < U ( J ; 1 ) -WNE<AN < 1 > > .♦ < W ( J y ^ ) -WNEAN <>!>)> /N 
C0W2W3=:C0W2W3+< (W( Jy2>-WMEAN<2) >t.(W< Jr3)-WNEAN<3> ) >/N 
C0W2W 'i =C0W2W ^ -f ( < U < J y 2 ) -WNEAN ( 2 > > * ( W < J y >1 ) -WNEAN < ^. > ) ) /N 
C0W3W^r;C0W3W^.+( <W( Jy2)-WNEAN<3) >S'-< W< Jr <1 )-WNEAN ( 'I > > >/N 
COE 1 B2-C0B 1 B2+ ( < E < J y i ) -ENEAN ( 1 ) ) « ( E < J r 2 > -BNEAN ( 2 > > > /N 
CGBlB3n:C0ElB3+( (E< Jy 1 > -BNEAN <1 > Jr3)-ENEAN <3 ' ' ' 

COB 1 B 'I ^COB 1 B 1 + ( ( B < J y 1 ) -BNE AN < 1 ) ) S: ( E ( J r >! > -ENE AN < ^ > > > /N 
CGB2B3nC0E2E3+( <E( Jr2)-BNEAN<2) ).'!c<E( Jr3>-ENEAN<3) > >/N 
C0B2B'!=C0E2E^+( (B( Jy2)-BNEAN(2) )J!c(E( Jr ^ ) -BNEAN <'! ) > )/N 
C0B3B^=C0B3B^+( {B( Jy3)-BNEAN(3) >S:(E( Jr > -BNEAN < 'S > > >/N 
COVl V2.-C0 VI V2+ ( ( V ( J y 1 ) - VNEAN ( 1 > ) S'. < V < J r 2 ) - VNEAN < 2 ) > ) /N 
C0VlV2-C0VlV3+< <y< Jr 1>-VNEAN(1) >*('.'( Jy3)-VNEAN(3) ) )/N 
COVl V^.-C0V1 V^+ ( < V < J r 1 ) -VNEAN ( 1 ) > S ( V < J r n ) -VNEAN ( ^ ) ) > /N 
C0V2 V3.-CGV2 V3+ < < V < J r 2 ) -VNEAN ( 2 ) ) S'. < V < J r 3 > - VNEAN ( 3 ) ) ) /N 
C0V2V'l.-CaV2V^+< <V< Jr2>-VNEAN<2> ).t.(V( Jr ^ >-VNEAN( ^ ) > >/N 
CO V3 V T --C0 V3 V 'J ■!■ < < V ( J r 3 ) - VNE AN ( 3 > ) S: ( V ( J r 'I ) - VNE AN ( ^ > ) > /N 
COR 1 R2.-CGR 1 R2+ < ( R C J r 1 > -RNE AN ( 1 > > S'. C R ( J r 2 > -RNE AN < 2 > ) ) /N 
COR 1 R3=C0R 1 R3+ < < R < J r 1 > -RNE AN < 1 > > S: < R ( J r 3 > -RNE AN ( 3 ) > ) /N 
COR 1 R ^.-CORl R -1 + < ( R < J r 1 > -RNE AN < 1 > ) S'. < R < J r 1 > -RNE AN ( ^ > ) ) /N 
C0R2R3-C:GR2R3+ < ( R ( J r 2 > -RNE AN (2>>S:(R<Jr3> -RNE AN ( 3 ) > ) /N 
C0R2R -1 .-C0R2R 1 + ( ( R < J r 2 ) -RNEAN < 2 ) ) S; { R ( J r A ) -RNEAN <'!)>) /N 
C0R3R'l--C0R3R^ + ( (R( Jr3>-RNEAN(3) ) S'. (R < J r -1 > -RNEAN < -1 > > >/N 
CONTINUE 

• 4 * -T**#* T* •»••»• •!•••• T'T* 

CALCULATION OF STIFFNESS NATRICES 


CALL WSTIFF < WVAR r WKNAT ) 

nAI I MAXTI.I / Ulc-MAT - . riCT ^ 

PAI I DCTTC'Cr / DIIAD . X>k*WAT \ 

PAI I WATTXI ^ r^k MAT . ^ « riCTT N 
HA! I I ICTTCrc: ^ I II lAD • I Ik'MAT \ 

PAI I WATTKI/ I IK^MAT. ^ ^r^CTT % 
nAI I OeTTCrer ^ C.| lAO - Dk-MAT \ 

W*t 4 -.k. I^W* A. . * 

rAf I MATTkl I' Dk' MAT - ^ .riCTT ^ 

tli.t a, » % % t \ $ %t mt i * # l#4U'k.. • ^ 

UpTTCr / 07 «TrO A \ 

•*vr%ai. i ita. WA> •a’ 

ernpM AT f 1 AV . ' Uk*MAT / % 

* ^ 

UDTTC/ 07 « 70 c:> //MkMAT/T. * 77 ** ? > r Isr.l - >» > 


pnOM AT / 1 AV « / OWMAT ^ ' 

« W« %«aj«4 A.*l%lfi»4 # 

WRI TE ( 23 y 325 ) ( ( BKNAT ( I r J) r J=1 y '5 > r I"! r 'J ) 
WRITE{23y327) 

pnow AT / 1 AV - ' I tk"M AT / \ 

WRITE < 23 r 325 ) < ( VKMAT C I r J ) r J«»l r ^ ) r 1*^1 r ) 
WRITE<23r32S> 



• '/ o o 

w 


c 

c 

c 


r 

Cl COO 
C 

r i AA 

^ w .fc 


00 

* t 

1 AA 

J. w v 


O A 

W w 


•7A 

t \f 


c 

c 

c 


11 


FORMAKIOX? 'RKMAT' > 

URITE(23y325) < < RKMrtT ( I ? J> J-i - ^ , t-j . » % 
F0RMriT(lXr'!(F15,6j3X)) 


REftDIK'G THE TEST SAMPLES AMD THEIR ANALYSIS 

Iv' — OAA ■ < • •• 

i % 4Ui vr w 

READ (21 yS') ( ( Y( I y J> 7 I«:l rK) ? J=?l 7 ^ ) 

UC'T TC / t AAA•^ 

FORMAT ( IX ; ' GBSERVATIDM SAMPLES ' ) 

UR I TE < 2 ? 1 OO 1 > < ( Y ( I E A ? I R A > 7 1 RA= 1 T -1 > 7 1 S A= 17 200) 

CnOMAT / 1 V . 7 T /» . TV % \ w w . 

I.IDTTC* /O/l « OO \ 

I# / / i 

crnoMAT i czv ^iia/'-/\v 

W D TTC / O A - i A A % 

IvU. •> # 

F ORMAT < 5X 7 ' SAMPLE MG ,' 73 X 7 ' CLASS NO , ' 7 1 OX 7 ' UU ' 7 1 2X 7 ' BU ' 7 
110 X 7 'L'U" 71 ^X 7 'RU' 7 IOX 7 'DECISION' ) 

UDTTET (T OA . OO \ 

•# # * 0 

DO 70 I^lrK 
II-I 

Tin OA l-i • A 

JL*w ww ^ 0 h, 0 t 

MU/ I\ 5 =.V/T. IV— UMCTAKI/ IV 
XfT/r IV~V/T« IV— OMCTAKI/ IV 

\ 0 k. 0 0 

nil/ I V —Y ^ T « IV— IIMITAVI/ IV 

w-<0%<».00‘^i%0u0%j0 

C'p/ IV=iY^T« l\— PWCTAXI/ IV 

£;n>viT JRiy£ 

call" PRODUT ( yK.MAT 7 UW I WU > 

nAi I oonni IX / ok’MAT - t>t> « oil V 
pAi I c»i:*nriiiT/iik*wAT«iiii«iiii V 

CALL PRODUT ( RKMAT 7 RR 7 RU ) 

TC 7 TT I C Tf'l ACO-i 

X M \ X X 1r •mmtm, ^ %0 0 <Ji.W4m.IIWW****. 

Tcr/TT OT c:AV Tri ACC-iO 

At \ X 0 k. r $ 9 %.0 %0 0 AWlMtlWW *~*<tM> 

TCT/TT pr iAAV TAI ACC-r* 3 : 

A4 \AAVWt * ,0. %0 90 0 AwLmI tw W** W 

IF<II.GT.150) ICLASSr^. 

CALL UM I MUM ( 1 1 7 1 CLASS 7 UU 7 BU 7 VU 7 RU 7 1 BES I > 

rnMTTI 1 .ll ic 

WWt%t At\W4«. 

UDTTir / O-A - OO V 

TTi««t.taM.%<t_ I# 0 0 0 

STOP 

END 

St/ ^t/ SI/ St/ st/ ^i/ SI/ Si/ St/ S/ S/ ^/ S^ si/ S^ ^/ S/ St/ St/ SI/ ^/ ^/ Si/ ^i/ ^ 

SUBROUTINE TO FORMAT A MATRIX UHOSE ELEMENTS ARE INVERSE 
OF VARIANCES AND COVARIANCES 

^ S/ S/ "S/ ^i/ St/ ^i/ SI/ Sf/ St^ ^1/ ^/ St/ Si/ ^X^ ^X* S/ ^X ^1/ Si/ S/ X St/ SI/ ^X 'X X' *X X* X^ *X X' ^X St^ 'X x^ x^ 

Cl it>c»ni iTTMcr ucTTcrc: /i.ii iac* - uk-mat v 

WW4U»4\WWtdkt%»M. WW*.4k« • <kVrvi*l%tW4\t»t4t« 

DIMEMS lOM WKMAT ( 7 '! } r yy<^,R ( 'I ) 

COMMON/AREAl /C0W1U2 7 C0U1U3 7 COUIU 7 C0U2W3 7 CGU2U'1 7 C0W3U 'i 
DO 11 I7=l7^ 

MK'M AT / T . T V «sl.II I AO /TV 

9 0 U 0 a. 0 99 9 t t* 9 1 .*.# 

CONTINUE 

WKMAT<l72)=C0UlU2 
NKMAT (271) =UKMAT < 1 7 2 ) 

UKMAT(l73)<=C0UlU3 
UKMAT ( 3 7 1 > =UKMAT (17 3) 



11 


11 


UKHAT(2r3)=C0W2W3 
WKM AT ( 3 ; 2 > •- WKMAT ( 2 ; 3 > 

U'KMAT ( 2 - 'I > -C0W2W1 
WKMAT < n » 2 ) -.WKMAT ( 2 ? 'S ) 

UKMAT<3r'5):cCOW3Wn 
WKMAT <'5 73) --WKMAT (37'!) 

RETURN 

END 

•!• •!• ««•»•* *|*««*»g«*(*«g* «4««g« •7* ■•• *i* •»•«*• *!• «$•*••••••■**§* •!•••'» *1* •«••••♦•* *|«*0**t^ •#•*••••• ►#••■• •!•*«• ••*•»• ••*»•» • » 

SUBROUTINE BSTIFF <ByAR7 EKMAT) 

nTMFMCTnM ck- MAT^ a .a \ . oiiAtj/ a y 

COMMGN/AREA2/COEiE27COBlB3 7 COBlE' 57 COB 2 E 3 7 CaE 2 B'!rCGB 3 B'! 

DO 11 I'=i7'l 
BKMAT(l7l):=ByAR(I) 

CONTINUE 

BKMAT(l72).-CDElB2 
EKMAT (271) -EKMAT (172) 

UK.'MAT / i . T \ — f-noi Vil 

BKMAT (371) -BKMAT (173) 

BKMAT(l7'))=C0BlB'! 

Cik-M AT / A - 1 ^ — trV'M AT 7 1 . A > 

*.*1 %t ll i I V I # ^ V* tt • 1 \ A ^ I # 

t>K‘*>wt AT / O - "2 \ 

BKMAT (372) :^KMAT (273) 

BKMAT(27'!)=C0B2E'! 

BKMAT ( A 7 2 ) --BKMAT ( 2 7 A ) 

EKMAT (37 '!)-.C0E3B'! 

BKMAT ('173) --EKMAT (37'!) 

RETURN 

END 

«|» »!• «l« *!• •I#*!* «g- «J**t«*t* •i**!**** *•* *1 ••• •#»*■•*••• 

SUEROUT I NE VST I FF ( WAR ? VKM AT > 

riTMTKICT nxi llk*MAT ( A •. A ) r ) 

COMMO^^AREA^'Cwiw 
BO 11 I-l?'! 

I Ik-M AT / T . T ^ r^ni lAD /TV 
rnKiT TKii ic 

wwt% ■ dk.t«w*>. 

llk* MAT/1 « O V —rni M I lO 

\ / ••ww T A V A. 

VKMAT ( 2 7 1 ) »VKMAT ( 1 7 2 ) 

1 Ik-'MAT 7 1,1 \ — oni 11 I IT 

k.d.#W / •'"•w W * * T w 

Ilk MAT/T^I V zr,« n-'h<AT ( 1 r 3 ) 

V41.I4III \ \ a.7W# 

I lk*M AT f ^ ^ A V —r ni !1 I \A 

Vt«l4llt » A W » A » « 

VKMAT ('!-!) -VKMAT ( 1 r '! ) 

VKMAT (273) =:C0 V2V3 
VKMAT (372 ) WKMAT (273) 

VKMAT(2 7'!).-C0V2V'! 

VKMAT ( 'i 7 2 ).-.VKM AT (27-1) 

VKMAT (37'!) =C0V3V'! 

VKMAT (-173) =VKMAT ( 3 7 'I ) 

RETURN 

END 

*/*********.*****.*-*.*.**.*.*.**.*.***^.^.^.*-'i'-’i^*^-'^-****^”**^**^'^'^‘^^^’^‘^' 

SUBROUTINE RSTIFF( RVAR 7 RKMAT ) 

DIMENSION RKMAT ( -1 7 -1 ) 7 RVAR ( '! ) 



60 


100 


c 

c 

c 

c 

c 

c 

Q 

c 


c 

10 

15 

20 

C 

30 

^0 

<5 

50 

wO 

70 

SO 

o>=: 

w w 

90 

95 

1 riA 

105 

110 

C 

1 -TA 

* W W 
1^0 
150 

1 ZA 

170 

JU_ V 

ox A 

W W 

0*7 A 

^ / vr 

c 

210 

TOA 

W Ax W 

“7*7A 

WW W 

*7-^ A 

W * 

2S0 

r 


r^^*TTKI» ic* 

URITEiS'Jyl 
FORMAT (SX 7 
112 Xt 12) 


00 > 1 1 7 1 CLASS 7 UU 7 BU 7 K'U 7 RU 7 IDESI 
I3 7 9X7l2 75X7F12.>)73X7F10,37 5X7F12,'‘ 


RETURN 


/iv.c-irr A . 

»/%y« •t.wv 


END 


SUBROUTINE FOR MATRIX INVERSION 

SUBROUTINE MATIN( A^NyDETERM) 
A»CO-EFFICIENT OF ORDER N 
E“=VECT0R OF ORDER N 


M:=IF M IS SET TO ZERO 7 ONLY INVERSEIS COMPUTED 
DETERM«VALUE OF DETERMENENT RETURNED 


DIMENSION A<^.7^) -IPIVOK^) yINDEX<^.r2) 

EOUI VALENCE < I ROW y JROU ) 7 ( ICOLUM 7 JCOLUM > - ( AMAX 7 T - SWAP > 

INITIALIZATION 

DETERM=1,0 

DO 20 J=l7N 


IPIV0T(J)«--0 

SEARCH FOR PIVOT ELEMENT 


DO 550 I-I 7 N 
AMAX»0,0 
DO 105 J«l7N 

IF(IPIVOT< J)-l >£07105760 
DO IOC K^lrN 

IF(IPIVOT<K)-l)S07l00 77'!0 

IF<AMAX-ABS{A<J7K)))S57lOOylOO 

IROW--J 

ICOLUM--=K 

AMAX=ABS<A( JyK) > 

CONTINUE 

CONTINUE 


I F I VOT (ICOLUM) *=I P I VOT ( ICOLUM ) + 1 

INTRECHANGE ROUS TO PUT PIVOT VECTORS ONDIAGONAL 
IF< IR0U-ICCLUM)1''.0 7 2607 1 TO 

riC-TcroM— —riETxroM 

I *..1 %« •**■ %1 t 

DO 200 L^lrN 

C'l.IAC-siA / TCni.l _ I \ 

w W i i j “7 \ / W / 

A < I RGU 7 L ) --A ( I COLUM 7 L ) 

A(IC0LUM7L)=^WAP 

TMnCV/T.I ^3;Tl^nU 

'•’NDEX (17 2) =1 COLUM 

DIVIDE PIVOT ROU BY PIVOR ELEMENT 
PI VOT- A ( I COLUM 7 1 COLUM ) 

DETERM=DETERM*PIVOT 
A < ICOLUM? ICOLUM)-! ,0 
DO 350 L=i7N 

A < I COLUM 7 L ) «A < I COLUM 7 L ) /P I VOT 
REDUCE NON PIVOT ROUS 



•70A 

w w V 

390 

'lOO 

“120 

■^,30 

^50 


r 

w 

c 


DO 550 L1:-1tN 

^.00^5507 'JOO 
T’^ritLl jIcoLUM) 
r>.LlrlCQHJM)r:0,0 
DO nSO L=1 


Kl 


0.<LlrL>=:A<LlTL>“/^(ICOLU«rL)ScT 

CON'TI^uic 


r 

w 

IMTERChamCE the COLUMNS 

500 

DO 710 1 = 1 , N 

510 

L-N+l-i 

^OA 

w*-. W 

IF<INr!E:v<i_7l)-INDEX(L-2 

Z*7A 

WW W 

JR0Wr-i.^,ccx<L7l) 

5^10 

JCOLUM~'iKiriEX(L» 2 ) 

550 

DO 705 K=lyN 

-t i A 

w V/ ^ 

SAP5=A<H:7 JROU) 

/."7A 

W # vr 

( K 7 JRow ) =A ( K 7 JCOLUM > 

”7 AA 

# vr w 

A<K7 JCQL'JM>-SWAP 

705 

CONTINUE 

710 

CONTINUE 

DO 11 

IF(IPIUot(K> ,NE.1)G0 TO 

11 

CONTINUE 

RETURN 

12 

URITE(227P91) 

001 

/ # J. 

FORMAT' /30X 7 'MATRIX IS 

-7ArS 

# 1 w 

RETURN 


E.ND 

^ Nl; J; ^ ^ ^ ^ ^ ^ 4 , 2^ ^ ^ 2*^ 2}; .'Ji ^ ]f5c ^ J 



